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ABSTRACT 

Absolute  proper  motions  are  determined  for  stars  and  galaxies  to  V  =  17.5  over  a  450  deg2  area  that  encloses  both 
Magellanic  Clouds.  The  proper  motions  are  based  on  photographic  and  CCD  observations  of  the  Yale/San  Juan 
Southern  Proper  Motion  program,  which  span  a  baseline  of  40  years.  Multiple,  local  relative  proper-motion  measures 
are  combined  in  an  overlap  solution  using  photometrically  selected  Galactic  disk  stars  to  define  a  global  relative 
system  that  is  then  transformed  to  absolute  using  external  galaxies  and  Hipparcos  stars  to  tie  into  the  ICRS.  The 
resulting  catalog  of  1 .4  million  objects  is  used  to  derive  the  mean  absolute  proper  motions  of  the  Large  Magellanic 
Cloud  (LMC)  and  the  Small  Magellanic  Cloud  (SMC);  Qju  cosS,  Ps)lmc  —  (1-89,  +0.39)  ±  (0.27, 0.27)  mas  yr~l 
and  0racos  5,  jtUsmc  =  (0.98,  - 1.01)  ±  (0.30, 0.29)  mas  yr“*.  These  mean  motions  are  based  on  best-measured 
samples  of  3822  LMC  stars  and  964  SMC  stars.  A  dominant  portion  (0.25  mas  yr_1)  of  the  formal  errors  is  due 
to  the  estimated  uncertainty  in  the  inertial  system  of  the  Hipparcos  Catalog  stars  used  to  anchor  the  bright  end 
of  our  proper  motion  measures.  A  more  precise  determination  can  be  made  for  the  proper  motion  of  the  SMC 
relative  to  the  LMC;  (MacosS.  A9)smc-i.mc  =  (—0.91,  —1.49)  ±  (0.16, 0.15)  mas  yr"1.  This  differential  value  is 
combined  with  measurements  of  the  proper  motion  of  the  LMC  taken  from  the  literature  to  produce  new  absolute 
proper-motion  determinations  for  the  SMC,  as  well  as  an  estimate  of  the  total  velocity  difference  of  the  two  clouds 
to  within  ±54  km  s_1 .  The  absolute  proper-motion  results  are  consistent  with  the  Clouds’  orbits  being  marginally 
bound  to  the  Milky  Way,  albeit  on  an  elongated  orbit.  The  inferred  relative  velocity  between  the  Clouds  places  them 
near  their  binding  energy  limit  and,  thus,  no  definitive  conclusion  can  be  made  as  to  whether  or  not  the  Clouds  are 
bound  to  one  another. 

Key  words;  astrometry  -  catalogs  -  Magellanic  Clouds  -  proper  motions 
Online-only  material:  color  figure 


1.  INTRODUCTION 

The  Magellanic  Clouds  have  provided  astronomers  with  a 
wide  variety  of  information,  from  the  small  to  the  large  hierarchy 
of  objects  in  the  universe.  They  are  the  first  step  in  the  cosmic- 
distance  ladder,  as  well  as  a  proxy  for  the  low-metal  gas-rich 
galaxies  assembled  in  the  early  universe,  they  have  zones  of 
recent  massive  star  formation,  young  1-3  Gyr  globular  clusters, 
and  a  variety  of  old  pulsating  stars  and  planetary  nebulae. 
They  contain  substantial  amounts  of  gas  and  dust  in  a  violent 
interstellar  medium,  harbor  the  closest  supernova  in  recent 
years,  and  are  being  used  as  a  testbed  for  dark  matter  searches. 

The  Magellanic  Clouds  are  also  the  prime  example  of  a 
galaxy-galaxy  interaction,  based  on  several  lines  of  evidence, 
close  enough  to  be  studied  in  detail;  an  apparently  starless 
Magellanic  Stream  trails  the  Clouds  (Mathewson  et  al.  1974; 
Penston  1982);  a  bridge  of  gas  and  stars  connects  the  Clouds 
(Hindman  et  al.  1963;  Irwin  et  al.  1985);  a  gaseous  Leading 
Arm  precedes  the  clouds  (Putman  et  al.  1998);  an  Hi  envelope 
surrounds  the  whole  system  and  a  collection  of  high-velocity 
clouds  seems  to  be  “raining”  over  the  Galaxy  (Olano  2004). 
Despite  such  convenient  observational  circumstances,  on  the 
theoretical  side,  no  dynamical  model  has  been  able  to  reproduce 
all  these  phenomena  simultaneously. 
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It  is  widely  accepted  that  the  complexity  and  intricacies  of 
the  Magellanic  Clouds’  external  and  internal  features,  have 
been  largely  determined  by  the  orbit  they  have  followed  in  the 
past  few  Gigayears.  Due  to  their  large  distance,  about  50  and 
60  kpc  to  the  Large  Magellanic  Cloud  (LMC)  and  the  Small 
Magellanic  Cloud  (SMC),  respectively,  only  radial  velocities 
were  precise  enough  to  provide  some  assessment  of  their  spatial 
velocity.  In  fact,  line-of-sight  measurements  of  the  Clouds  began 
about  a  hundred  years  ago  (Wilson  1935),  while  proper-motion 
measurements  of  a  useful  accuracy  were  possible  only  in  the 
1990s. 

The  first  proper-motion  results  (Jones  et  al.  1989;  Tucholke 
&  Hiesgen  1991;  Bastian  et  al.  1993;  Lin  1993;  Kroupa  et  al. 
1994;  Jones  et  al.  1994;  Irwin  et  al.  1996;  Kroupa  &  Bastian 
1997a,  1997b,  1997c;  Anguital999)  based  on  plate  and/or  CCD 
data,  were  compatible  with  a  picture  in  which  the  Magellanic 
Clouds  were  bound  to  each  other  and  to  the  Milky  Way.  Such 
scenario  relied  heavily  on  the  fact  that  the  Galactic  gravitational 
potential  used  (isothermal  sphere)  yields  such  results  by  default, 
and  proper-motion  errors  were  not  small  enough  to  refine  the 
tangential  velocities. 

In  the  past  10  years  though,  investigations  yielded  quite  a  va¬ 
riety  of  results  (Anguita  et  al.  2000;  Drake  et  al.  2001;  Pedreros 
et  al.  2002, 2006a,  2006b;  Momany  &  Zaggia  2005;  Kallivayalil 
et  al.  2005,  2006a,  2006b,  2009;  Pedreros  2006;  M6ndez  et  al. 
2006;  Piatek  et  al.  2008;  Costa  et  al.  2009).  Some  (Anguita  et  al. 
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2000;  Momany  &  Zaggia  2005)  were  found  to  have  unknown  or 
important  systematic  errors.  More  interestingly,  Hubble  Space 
Telescope  (//S7)-based  results  (Kallivayalil  et  al.  2005,  2006a, 
2006b,  2009;  Piatek  et  al.  2008)  coupled  to  more  modem  and 
cosmologically  inspired  dark  matter  Halo  models,  suggest  that 
the  Clouds  are  not  bound  to  the  Galaxy,  opposite  to  the  long-held 
paradigm. 

Twenty  years  have  passed  since  the  very  first  proper-motion 
measurements  of  the  Clouds,  and  it  is  only  now  with  the  Yale- 
San  Juan  Southern  Proper  Motion  (SPM)  program — briefly 
explained  in  Section  2 — that  for  the  first  time  a  wide-field 
astrometric  proper  motion  survey  of  the  Magellanic  System  is 
finally  completed.  All  the  first-epoch  (early  1970s)  and  part  of 
the  second-epoch  (early  1990s)  SPM  material  used  in  this  work 
are  photographic  plates.  Their  processing  is  briefly  summarized 
in  Section  3,  but  a  more  detailed  explanation  can  be  found  in 
T.  M.  Girard  et  al.  (2010,  in  preparation).7  A  substantial  part 
of  our  second-epoch  data  comes  from  SPM  CCD  observations. 
A  short  explanation  of  the  data  acquisition,  quality  control,  and 
processing  is  explained  in  Section  4.  To  supplement  our  second- 
epoch  plate  data  with  CCD-quality  positions,  we  have  included 
mean  positions  at  Julian  epoch  2000.0  from  the  UCAC2  catalog 
(Zacharias  et  al.  2004).  Section  5  explains  how  these  data  were 
selected  and  included  in  this  investigation. 

In  contrast  to  other  SPM  reductions,  as  explained  in  Section  6, 
relative  proper  motions  measured  in  CCD-size  fields  of  view 
were  combined  into  a  single  common  extended  and  accurately 
defined  global  system.  Although  our  zero  point  accuracy,  i.e., 
how  well  our  reference  frame  is  linked  to  the  International 
Celestial  Reference  System  (ICRS),  is  ultimately  limited  by 
Hipparcos  accuracy,  our  very  precise  relative  proper  motions 
over  the  whole  field  of  view  (FOV),  enabled  us  to  measure  the 
proper  motion  of  the  SMC  with  respect  to  the  LMC,  at  a  pre¬ 
cision  comparable  to  the  quoted  errors  of  space-based  proper 
motions. 

It  is  this  capability  that  we  exploit  to  obtain  new  measure¬ 
ments  of  the  proper  motion  of  the  SMC  based  on  previously 
published  LMC  proper  motions.  Section  7  contains  the  main 
results  of  this  paper  regarding  the  proper  motion  of  the  Clouds, 
absolute  and  relative.  Section  8  has  a  discussion  of  the  implica¬ 
tions  of  our  results  on  the  current  understanding  of  the  dynamics 
of  the  Magellanic  System.  Finally,  Section  9  states  the  conclu¬ 
sions  of  this  paper,  and  future  plans  already  in  consideration  to 
improve  the  current  results. 

2.  THE  SPM  PROGRAM 

This  investigation  is  part  of  the  SPM  program,  a  joint  venture 
of  Yale  University  and  Universidad  Nacional  de  San  Juan  in 
Argentina.  The  SPM  program  was  initiated  in  the  early  1960s 
by  D.  Brower  and  J.  Schilt  as  a  joint  enterprise  of  the  Yale  and 
Columbia  Universities  (Wesselink  1 974).  The  goal  of  the  SPM 
program  is  to  provide  absolute  proper  motions,  positions,  and 
BV  photometry  for  the  Southern  sky  to  a  limiting  magnitude  of 
V  ~  18. 

The  SPM  program  makes  use  of  the  Yale  Southern  Observa¬ 
tory’s  double  astrograph  at  Cesco  Observatory  in  the  foothills  of 
the  Andes  mountains  in  El  Leoncito,  Argentina.  This  telescope 
consists  of  two  51  cm  refractors,  designed  for  photography  in 
the  blue  and  yellow  bands,  respectively.  The  first-epoch  survey, 
taken  between  1965  and  1974  was  made  on  glass  photographic 
plates,  exposed  simultaneously  in  blue-yellow  pairs  and  always 

7  See  also  http://www.astro.ya1e.edu/astrom/spni4cat/spm4.html 


centered  on  the  meridian.  The  plates’  FOV  extends  over  an  area 
of  6°.  3  x  613.  The  sky  south  of  S  —  —17°  was  observed  in  the 
first-epoch  period. 

Second-epoch  SPM  plate  observations  were  begun  in  1988. 
By  the  mid  1990’s,  with  only  a  third  of  the  second-epoch 
survey  completed,  Kodak  discontinued  the  production  of  the 
photographic  plates.  In  2000,  with  funding  from  the  NSF,  a 
CCD  camera  system  was  installed  on  the  double  astrograph  to 
replace  the  photographic  plate  holders.  A  PixelVision  4K  x 
4K  CCD  camera  (0?94  x  0194  FOV)  was  placed  in  the  yellow 
lens  focal  plane,  and  an  Apogee  Ap-8  IK  x  IK  (0137  x  0137 
FOV)  was  fitted  in  the  blue  focal  plane.  In  2004,  the  Apogee 
IK  camera  was  upgraded  to  an  Apogee  Alta  E42  2K  x  2K 
(0142  x  0142  FOV)  with  funds  from  the  Argentine  CONICET. 
In  order  to  achieve  a  limiting  magnitude  similar  to  the  first-epoch 
plate  material,  the  CCD  survey  consists  of  2  minute  exposures 
in  both  cameras. 

In  the  past  decade,  catalogs  of  the  SPM  program  covering 
various  parts  of  the  sky  have  been  published,  as  the  second- 
epoch  material  became  available  for  its  astrometric  reduction. 
Catalogs  SPM1  (Platais  et  al.  1998),  SPM2  (van  Altena  et  al. 
1 999),  and  SPM3  (Girard  et  al.  2004)  are  based  on  photographic 
plates  only,  in  both  first  and  second  epochs.  The  plates  were 
scanned  either  with  the  Yale  PDS  (using  input  lists  of  selected 
objects)  or  the  USNO  Precision  Measuring  Machine  (PMM)  for 
the  whole  plate.  Also,  two  different  centering  algorithms  have 
been  used  to  measure  the  image  centers  in  the  scans,  the  Yale 
two-dimensional  Gaussian  fit  from  Lee  &  van  Altena  (1983),  or 
the  USNO  circular  fit  from  Monet  et  al.  (2003). 

Since  2004,  regular  CCD  observations  have  been  carried  out 
to  finish  the  second-epoch  survey  of  the  SPM  program.  By  2008 
December,  the  survey  was  effectively  completed  for  the  sky 
south  of  5  =  —20°.  Subsequently,  the  SPM4  catalog,  based  on 
all  available  plate  and  CCD  data,  was  completed  in  late  2009  and 
is  currently  available,  (T.  M.  Girard  et  al.  2010,  in  preparation). 
SPM4  includes  ~100  million  objects  south  of  S  =  —20°, 
brighter  than  V  ~  18.  Many  of  the  data-processing  procedures, 
software,  and  protocols  developed  in  this  investigation  were  also 
used  in  the  construction  of  the  SPM4  catalog. 

3.  THE  PLATE  DATA 
3.1.  Observations 

The  SPM  survey  fields  are  on  5°  centers  in  declination  and  a 
maximum  of  5°  separation  in  right  ascension,  providing  at  least 
a  full  degree  of  overlap  between  adjacent  plates.  Each  17  inch 
x  17  inch  plate  covers  an  area  of  6? 3  x  6? 3  (55"  1  mm-1  plate 
scale)  and  consists  of  a  2  hr  and  a  2  minute  offset  exposure.  All 
observations  were  made  with  a  wire  grating  over  the  objectives, 
producing  measurable  diffraction  images  out  to  third  order. 
The  grating  constant  is  3.8  mag,  thus,  along  with  the  offset 
short  exposure,  effectively  increasing  the  dynamic  range  of 
each  plate  allowing  measurement  of  external  galaxies  and  bright 
Hipparcos  stars.  Fields  were  observed  simultaneously  in  blue 
and  yellow  passbands,  and  there  are  some  fields  with  repeated 
blue  and/or  yellow  plates  from  the  same  epoch.  See  Girard 
et  al.  (2004)  for  a  more  complete  description  of  the  SPM  plate 
material. 

Table  1  lists  the  SPM  plates  used  in  this  investigation.  During 
the  course  of  this  research,  it  was  found  that  some  plates  yielded 
unusually  deviant  results.  A  visual  examination  of  the  suspect 
plates  revealed  that  the  stellar  images  suffered  from  significant 
defects,  possibly  caused  by  poor  guiding,  polar  misalignment, 
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Table  1 

SPM  Plates  Used  in  This  Paper 


SPM  Field 

(a.  5 ) i 950 

SPM  Plates 

1st  Epoch 

2nd  Epoch 

Plate  Number 

Date 

Plate  Number 

Date 

028 

(00h00m,  -75°) 

0829B/Y 

03/NOV/72.84 

1354/BY 

01 /NOV/94.84 

029 

(01h00m,  -75°) 

0834B/Y 

06/NOV/72.85 

030 

(02h00m.  -75°) 

0826B/Y 

01 /NOV/72.84 

031 

(03h00m,  -75°) 

0830B/Y 

03/NOV/72.84 

032 

(04h00m,  -75°) 

0950B/Y 

26/OCT/73.82 

033 

(05h00m,  -75°) 

0851B/Y 

06/JAN/73.02 

1357/BY 

30/NOV/94.92 

1373/BY 

20/OCT/95.80 

034 

(OS^OO111.  -75°) 

0860B/Y 

25/JAN/73.07 

052 

(OO^OO"1,  -70°) 

0821B/Y 

06/OCT/72.77 

1355/BY 

02/NOV/94.84 

053 

(00k48m,  -70°) 

0750B/Y 

24/AUG/71.65 

1371/BY 

23/SEP/95.73 

054 

(01h36m,  -70°) 

0751B/Y 

25/ AUG/7 1.65 

1372/BY 

20/OCT/95.80 

055 

(02k24m,  -70°) 

0969B/Y 

26/NOV/73.90 

056 

(03h12rn,  —70°) 

0815B/Y 

16/SEP/72.71 

057 

(04k00m.  -70°) 

0946B/Y 

06/OCT/73.76 

058 

(04h48m,  —70°) 

0S46B/Y 

01/JAN/73.00 

059 

(05h36m,  -70°) 

0854B/Y 

09/JAN/73.02 

1350/BY 

17/NOV/93.88 

060 

(06h24m,  -70°) 

0842B/Y 

13/DEC/72.95 

061 

(07h12ra,  —70°) 

0423B/Y 

20/FEB/69.14 

088 

(04h00ra,  -65°) 

0824B/Y 

15/OCT/72.79 

089 

(04h40™,  -65°) 

0835B/Y 

06/NOV/72.85 

090 

(05h20ra,  -65°) 

0853B/Y 

07/JAN/73.02 

1394/BY 

16/NOV/93.88 

091 

(Ob^O"1,  -65°) 

078  IB/ Y 

13/JAN/72.04 

0984B/Y 

27/JAN/74.07 

092  ' 

(06h40m,  -65°) 

0986B/Y 

28/JAN/74.08 

Notes.  Plates  O750B,  0751B,  1371B,  13S7Y,  and  1373Y,  were  subsequently  discarded,  as  visual  examination  of  them  revealed 
significant  deficiencies  that  would  negatively  impact  the  astrometric  reduction  at  their  locations. 


Figure  1.  SPM  fields  in  the  Magellanic  Cloud  area.  Gray  outlined  regions 
correspond  to  second-epoch  SPM  plate  observations.  The  small  squares  indicate 
second-epoch  CCD  observations.  The  black  outlined  area  is  the  VMC  region 
(see  Section  4.2).  A  subset  of  stars  from  the  USNO-A2.0  Catalog  is  plotted  in 
the  background  to  indicate  the  positions  of  the  LMC  (on  the  left)  and  SMC 
(on  the  right).  SPM  fields  are  on  nominal  5°  centers,  not  including  the  VMC 
field. 

or  poor  focus.  These  plates,  0750B,  0751B,  1371B,  1357Y,and 
1373Y,  were  therefore  discarded  for  the  research  presented  here. 
Coverage  in  these  areas  was  not  affected,  since  only  one  of  the 
two  plates  at  a  given  epoch  per  field  affected  was  discarded. 
Figure  1  shows  the  distribution  on  the  sky  of  the  plates  used 
in  this  work.  Twenty-two  (22)  SPM  regions  were  studied,  of 
which  seven  (7)  have  second-epoch  plates.  For  the  areas  with 
second-epoch  plate  data,  no  second-epoch  CCD  observations 
were  made. 


3.2.  Astrometric  Reduction  and  Photometric  Calibration 

The  SPM  plates  used  in  this  investigation  were  scanned  with 
the  PMM  of  the  USNO  Flagstaff  Station.  For  more  details  about 
the  PMM  setup  and  operations,  see  Monet  et  al.  (2003).  In  a 
collaboration  between  the  USNO-Washington  and  the  Yale  As¬ 
trometry  Group,  the  StarScan  reduction  pipeline  (Zacharias  et  al. 
2008)  was  modified  to  analyze  the  PMM  pixel  data  of  the  SPM 
plates  to  produce  a  list  of  detections,  image  centers  and  pho¬ 
tometric  indices.  The  astrometric  and  photometric  reductions 
then  proceeded  as  follows:  (1)  cross-identification  of  detections 
to  an  input  catalog,  including  (2)  recognition  and  identification 
of  central  and  higher  grating  orders,  (3)  photometric  calibration 
to  obtain  B  V,  (4)  correction  for  atmospheric  refraction,  (5)  cor¬ 
rection  for  Magnitude  Equation,  which  also  combines  grating- 
order  systems,  (6)  transformation  of  short-exposure  positions 
into  the  long-exposure  system,  and  (7)  astrometric  solution  into 
Tycho-2  to  obtain  (a,  S). 

The  input  catalog  referenced  in  step  (1)  is  a  compilation  of 
the  following  external  catalogs:  Hipparcos  (Perryman  &  ESA 
1997),  Tycho-2  (H0g  et  al.  2000),  UCAC2  (Zacharias  et  al. 
2004),  2MASS  point-source  and  extended-source  (Skrutskie 
et  al.  2006),  LEDA  galaxies  with  DENIS  measurements 
(Paturel  et  al.  2005),  and  the  QSO  catalog  of  Veron-Cetty  & 
Veron  (2006).  In  order  for  an  object  to  be  included  in  our  study, 
it  must  appear  in  one  or  more  of  these  listed  catalogs.  A  thorough 
explanation  of  all  these  procedures  can  be  found  in  T.  M.  Girard 
et  al.  (2010,  in  preparation).  After  the  above  processing,  one 
has  positions  (*,  y)  properly  calibrated  into  a  common  system 
within  each  plate,  with  computed  positional  errors  and  astro¬ 
nomical  coordinates  (a,  S )  on  the  ICRS,  as  realized  by  Tycho-2. 
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Figure  2.  Magnitude  distribution  and,  thus,  completeness  of  the  SPM  material 
evaluated  in  a  CCD-size  field  in  SPM  field  028  that  has  SPM  plate  data  (dark 
gray),  SPM  CCD  data  (light  gray),  and  UCAC2  data  (stitch  under  curve),  in 
comparison  with  the  2MASS  detections  (black)  in  the  field. 

3.3.  Evaluation  of  the  Plate  Data 

Well  measured  stars  on  the  plates  in  all  of  the  orders 
have  positional  errors  between  0.9 /xm  and  1.6 /xm  (50  mas 
and  90  mas).  These  errors,  which  only  assess  measurement 
uncertainties,  are  consistent  with  the  precision  expected  for  a 
good  centering  procedure,  based  on  previous  experience  with 
the  SPM  plates. 

A  single  final  position  ( x ,  y)  per  star  per  plate  is  obtained, 
from  the  positional-error-weighted-average  of  the  available 
measurements.  As  expected,  the  error  of  the  final  position 
and  magnitude  will  depend  on  the  number  of  grating  orders 
contributing  to  their  calculation.  If  average  measurements  from 
different  plates  were  later  averaged  to  obtain  a  final  number,  then 
other  errors,  random  and  systematic,  would  come  into  play,  and 
the  error  budget  of  this  final  result  should  include  the  additional 
sources  of  uncertainty. 

As  with  earlier  SPM  catalogs,  an  approximate  estimate  of  the 
relative  completeness  between  the  plate  data  and  Two  Micron 
All  Sky  Survey  (2MASS),  can  be  made.  Figure  2  shows  the  V 
magnitude  distribution  of  all  stars  detected  on  a  yellow  plate, 
compared  to  the  VJK  magnitude  distribution  of  all  2MASS  point 
source  stars  in  the  same  field.  Vjk  —  J  +  2.79(7  —  K)  is  an 
approximate  empirical  transformation  from  2MASS  JHK  to  V 
determined  by  Girard  et  al.  (2004)  and  found  to  be  valid  for  a 
relatively  wide  range  of  spectral  types.  It  can  be  seen  that  the 
SPM  plates  have  a  completeness  similar  to  that  of  2MASS  up 
to  V  =  17.5  and  a  falloff  after  that. 

In  general,  compared  to  previous  SPM  processing,  these 
plate  data  have  significantly  fewer  false  detections,  and  a  better 
correction  of  systematics  in  the  detected  positions  associated 
with  the  scanning  process. 

3.4.  SPM  Second-epoch  Plate  Data 

As  explained  in  Subsection  3.1,  seven  SPM  fields  have 
second-epoch  plates.  For  this  reason,  no  second-epoch  CCD 
observations  were  made  at  these  locations.  In  order  to  facilitate 
the  managing  of  files  and  software,  the  second-epoch  plates 
were  divided  in  CCD-size  fields,  to  emulate  the  overlap  scheme 
of  the  second-epoch  CCD  frames  (see  Section  4.2).  This  way,  all 
the  second-epoch  material,  regardless  of  its  type,  would  have  a 
uniform  format  and  structure,  for  programming  purposes.  From 


the  about  90  CCD  pointings  that  usually  cover  one  SPM  field, 
half  were  used  to  divide  the  blue  plate  and  the  other  half  to 
divide  the  yellow  plate,  in  such  way  that  these  yellow  and  blue 
fields  overlap  in  a  similar  way  as  do  the  real  CCD  frames. 

4.  THE  CCD  DATA 

4.1.  CCD  Cameras 

The  CCD  system  consists  of  four  cameras,  one  each  for 
the  blue-  and  yellow-optimized  lenses  and  two  focus-sensor 
cameras,  again,  one  for  each  lens.  The  CCD  camera  for  the 
Yellow  telescope  is  a  PixeiVision  camera  with  a  Cryotiger 
Chiller  cooled  (—85  °C)  4K  x  4K  Loral  chip  with  15  pm  pixels, 
which  translates  into  a  pixel  size  of  (1/83  and  a  total  area  of 
0/94  x  0?94  degrees.  The  pixel  size  is  well  matched  to  the  YSO 
site  where  the  seeing  conditions  usually  yield  an  image  FWHM 
of  2-3",  corresponding  to  3-4  pixels  per  FWHM.  This  is  optimal 
sampling  for  the  derivation  of  astrometric  image  centers,  based 
on  SPM  experience  with  digital  image  centering.  The  unthinned 
and  front-illuminated  Loral  chip  is  fitted  with  a  fixed  custom 
scientific  V-band  filter. 

The  Blue  telescope  was  first  fitted  with  an  Apogee  AP-8 
camera  that  utilized  a  thermoelectrieally  cooled  (— 40°C)  and 
back-illuminated  IK  x  IK  Site  chip  with  24  /xm  pixels,  which 
translated  into  a  pixel  size  of  1"32  and  a  total  area  of  22/57  x 
22/57.  In  2005  May,  an  upgrade  was  made  by  replacing  this 
camera  with  a  new  Apogee  Alta  E42  back-illuminated  2  K  x 
2  K  chip,  with  13.5  pm  pixels,  which  correspond  to  a  pixel  size 
of  0/'74  and  a  total  area  of  24/8  x  24/8.  Centered  on  the  same 
field  as  the  larger  PV  yellow  camera,  the  purpose  of  the  blue 
CCD  camera  is  to  provide  fi-band  CCD  photometry  for  the  stars 
that  fall  into  its  FOV. 

Data  from  the  yellow  PV  camera  were  used  for  both  astrom¬ 
etry  and  photometry,  while  the  blue  Apogee  and  Alta  cameras 
observations  were  used  only  for  the  photometric  calibration  of 
the  blue  plates. 

4.2.  Observations 

As  a  norm  for  the  SPM,  CCD  observations  are  done  always 
within  lh30m  of  the  meridian,  in  2  minute  exposures,  with  the 
wire  grating  placed  so  that  the  diffraction  pattern  is  at  about  45° 
from  the  E-W  line.  Normally,  an  E-W  orientation  is  ideal  to 
avoid  differential  color  refraction  effects  within  the  diffraction 
pattern,  but  in  this  case,  a  diagonal  orientation  prevents  the 
saturated  central-order  image  of  a  bright  star  from  spoiling  its 
grating  images  by  either  row  or  column  bleeding.  The  CCD 
pointings  conform  to  a  two-fold  overlap  coverage  scheme  for 
the  PV  camera,  as  shown  in  Figure  1 .  About  90  PV  CCD  frames 
cover  one  single  SPM  field.  Only  targets  that  did  not  have 
second-epoch  plate  data  were  observed  with  the  CCD,  except 
for  a  few  special  targets. 

A  total  of  1310  CCD  pointings  were  observed  for  this 
investigation.  Each  pointing  was  planned  to  be  observed  only 
once,  except  for  90  of  them  extending  over  a  6°  x  6°  field  around 
(a,  S)  =  (3h44m33s,  — 71°40'18"),  within  the  area  delimited  by 
the  bold  black  line  in  Figure  1.  This  area  corresponds  to  the 
non-SPM  field  for  the  variable  stars  in  the  Magellanic  Clouds 
(VMC)  study  done  from  1965  to  1968  by  A.  J.  Wesselink,  and 
comprises  seventy  60  minute  exposure  blue  plates  without  the 
wire  objective  grating,  reaching  a  limiting  magnitude  of  about 
18.  Numerous  repeated  CCD  observations  were  performed  on 
this  area  to  provide  suitable  second-epoch  observations  for  this 
material.  The  VMC  plates  were  not  used  for  this  work  because 
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they  require  their  own  special  reduction,  different  from  the  one 
used  in  normal  SPM  plates,  but  the  repeated  CCD  observations 
were  indeed  used  for  this  investigation.  The  VMC  CCD  targets 
have  at  least  13  good  observations  each,  with  some  of  them 
having  up  to  20  good  observations. 

The  criteria  used  to  qualify  an  observed  frame  as  acceptable 
are:  FWHM  ^  3"5,  limiting  magnitude  Vjjm  >  17,  and  standard 
error  <120  mas  from  an  astrometric  solution  into  the  UCAC2. 
If  a  frame  fails  any  of  these  limits  then  it  was  taken  again 
until  it  passed  all  of  them.  Nonetheless,  all  frames  regardless 
of  quality  are  saved  and  processed,  and  only  later  in  the 
astrometric  reduction  are  discarded,  if  they  prove  to  be  too  bad 
for  any  use.  In  total,  5422  CCD  frames  were  processed  for  this 
investigation. 


4.3.  Astrometric  Reduction  and  Photometric  Calibration 

Data  from  all  the  CCD  cameras  went  through  the  usual 
processing  to  calibrate  the  flux  detected  by  the  electronics,  for 
the  zero  charge  of  the  chip  (bias),  accumulated  signal  from  the 
electronics  dark  current  (dark)  and  different  response  to  light 
from  each  pixel  (flat).  Details  of  these  procedures  can  be  found 
in  T.  M.  Girard  et  al.  (2010,  in  preparation). 

Image  detection  on  the  processed  CCD  frames  is  done  using 
SExtractor  Version  2.4.4  (Bertin  &  Amouts  1 996),  from  which 
a  preliminary  centroid  and  an  aperture  instrumental  magnitude 
are  read.  SExtractor  centroids  are  then  used  as  input  positions 
to  compute  more  precise  centers,  based  on  the  Yale  two- 
dimensional  Gaussian  centering  algorithm  (Lee  &  van  Altena 
1983).  Typically,  more  than  90%  of  SExtractor  detections  are 
centered.  A  significant  reduction  in  the  rms  of  the  positions  for 
repeated  observations  of  the  PV  frames  was  found  when  using 
the  Yale-based  centers  as  compared  to  the  SExtractor-based 
centroids. 

The  process  of  transforming  the  Yale-based  centers  (x,  y)  and 
SExtractor-based  instrumental' magnitudes  minst  into  calibrated 
(a,  5)  and  BV,  for  the  CCD  frames,  follows  a  procedure  similar 
to  that  of  the  plates  (Section  3.2),  with  the  exception  that  the 
external  catalog  used  for  the  astrometric  reduction  of  the  CCDs 
is  UCAC2  (T.  M.  Girard  et  al.  2010,  in  preparation).  As  a  result 
of  these  procedures,  all  the  detected  positions  (x ,  y)  are  properly 
calibrated  within  each  CCD  frame,  with  computed  positional  er¬ 
rors  and  equatorial  coordinates  (a,  5)  in  the  ICRS,  as  realized  by 
UCAC2. 


4.4.  Evaluation  of  the  CCD  Data 

The  single-image  centering  precision  for  well  measured  stars 
(V  15)  in  the  CCD  frames  is  0.5  /rm  (25  mas),  worsening 
for  the  faintest  stars  where  it  reaches  about  2/zm  (100  mas). 
A  single  final  position  (x,  y)  and  magnitude  per  star  per  CCD 
is  obtained,  from  the  positional-enror-weighted-average  of  the 
available  measurements.  The  final  positional  and  photometric 
errors  will  depend  on  the  various  image  orders  contributing  to 
the  final  value. 

At  this  point,  each  star  in  each  CCD  frame  has  a  master  ID 
identification,  a  position  (x,  y)  with  errors,  a  calibrated  B  or 
V  magnitude  and  a  UCAC2-based  (a,  <$).  Because  stars  were 
identified  in  the  CCD  frames  in  the  same  way  as  in  the  plates, 
similar  charateristics  regarding  completeness  were  expected 
when  compared  to  2MASS  (see  Figure  2).  The  CCD  data  show 
in  general  a  completeness  magnitude  of  about  V  =  18. 


5.  UCAC2  CCD-POSITIONS  AS  SUPPLEMENT 
FOR  THE  2ND-EPOCH  PLATE  DATA 

As  seen  in  Figure  1,  some  fields  on  and  adjacent  to  the  Mag¬ 
ellanic  Clouds  have  only  photographic  plates  as  second-epoch 
material.  Given  the  lower  quality  of  plate  images  compared 
to  CCD  observations,  the  measured  proper  motions  in  these 
plate-only  fields,  will  have  significantly  larger  errors  than  those 
coming  from  the  combination  of  plate  and  CCD  data. 

In  an  attempt  to  counter  this  and  achieve  a  more  homogeneous 
quality  in  the  final  measurements,  we  decided  to  supplement 
the  second-epoch  plate  measures  with  epoch  2000  positions 
from  the  UCAC2  Catalog  (Zacharias  et  al.  2004).  These  are 
mostly  based  on  CCD  observations8  with  the  USNO  8  inch 
(0.2  m)  Twin  Astrograph  from  Cerro  Tololo  International 
Observatory  in  Chile.  The  UCAC2  data  were  collected  in  such 
way  that  it  mimics  our  second-epoch  fields.  The  precision 
of  the  positions  are  15  to  70  mas,  depending  on  magnitude, 
with  claimed  estimated  systematic  errors  of  10  mas  or  below. 
UCAC2  provides  only  crude  magnitudes  in  a  single  nonstandard 
bandpass  /?ucac2  between  V  and  R,  and  its  limiting  magnitude 
is  about  R  «  16. 

A  significant  number  of  faint  stars  in  our  second-epoch  SPM 
plates  do  not  have  a  counterpart  in  the  UCAC2  catalog.  This 
means  we  are  only  sampling  stars  about  V  <  16  for  these  SPM 
fields  when  using  UCAC2  positions.  UCAC2  completeness 
compared  to  2MASS  and  SPM  data,  can  be  seen  in  Figure  2. 

6.  OBTAINING  THE  PROPER  MOTIONS 

Given  that  a  substantial  part  of  the  second-epoch  SPM 
material  used  in  this  work  comes  from  CCD  frames,  previous 
SPM  procedures  for  obtaining  proper  motions  (used  on  first-  and 
second-epoch  plates),  could  not  be  straightforwardly  applied. 
Since  the  CCD’s  FOV  is  about  40  times  smaller  than  the 
plate’s  FOV,  the  number  of  reference  stars  available  to  measure 
relative  proper  motions  in  each  field  is  proportionally  smaller. 
A  simple  cut  in  magnitude,  as  in  past  SPM  reductions,  would 
result  in  too  few  reference  objects  per  CCD  frame;  in  particular, 
the  extragalactic  objects  needed  to  transform  relative  proper 
motions  into  absolute  ones. 

However,  if  we  select  reference  stars  belonging  to  some 
specific  population  of  the  Galaxy,  it  is  reasonable  to  assume 
that  they  have  a  mean  absolute  motion  along  an  extended  area 
on  the  sky9  that  can  be  parameterized  as  a  smooth  function 
of  (a,  5).  Moreover,  within  a  CCD  FOV,  their  mean  motion 
has  a  very  small  gradient,  if  any.  It  then  becomes  a  matter  of 
precisely  quantifying,  over  the  whole  FOV,  the  measured  mean 
relative  proper  motion  of  all  known  extragalactic  objects,  which 
is  simply  the  reflex  of  the  mean  absolute  proper  motion  of 
the  reference  stars.  Applying  such  a  function  to  the  measured 
relative  proper  motions  converts  them  to  absolute. 

The  layout  of  the  CCD  fields,  where  substantial  overlap  exists 
(20%  to  50%),  contains  a  wealth  of  information  that  can  provide 
linkage  of  the  reference  system  across  the  observed  region  of  the 
Clouds  that  is  limited  only  by  measurement  errors.  A  key  point 


8  Strictly  speaking,  the  supplemental  UCAC2  epoch-2000  positions  we  use 
include  up  to  a  few  years’  worth  of  UCAC2  proper  motions.  For  the  areas  of  . 
sky  in  this  study,  the  UCAC2  CCD  observations  were  made  around  1998.  Thus 
the  positions  we  employ  as  supplements  are  practically  those  of  the  UCAC 
CCD  program. 

9  Stars  behave  like  a  collisionless  system,  all  moving  under  the  influence  of 
the  same  general  background  gravitational  field,  therefore  we  can  expect  them 
to  have  a  global  smooth  distribution  in  their  velocities,  with  some  scatter 
around  a  mean  value  at  any  given  location. 
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Figure  3.  Sky  map  of  the  relative  proper-motion  reference  stars.  The  gray  curves  enclose  the  so-called  contaminated  areas. 


of  this  investigation  is  therefore  to  find  a  procedure  to  utilize 
this  large  overlap  to  produce  a  precise  global  reference  system. 

6.1.  Relative  Proper  Motions  in  the  Small  CCD  FOVs 

After  all  the  plate  and  CCD  data  processing  described  above, 
the  following  data  are  ready  for  the  measurement  of  proper 
motions. 

1.  First-epoch  positions  on  the  plates,  which  have  been  cor¬ 
rected  for  systematics,  as  much  as  Tycho-2  precision  and 
number  of  stars  available  allows. 

2.  Second-epoch  positions  on  SPM  CCD  frames,  SPM  plates, 
or  from  UCAC2.  CCD  positions  are  mostly  free  of  system¬ 
atic  errors.  Plate  positions  have  been  corrected  for  system¬ 
atics,  as  much  as  Tycho-2  stars  precision  and  number  of 
stars  allows. 

Henceforth,  we  will  refer  to  each  of  the  second-epoch 
CCD-size  frames  as  a  brick,  regardless  of  the  source  of  its 
data  (SPM  plate,  SPM  CCD  or  UCAC2).  For  each  second- 
epoch  brick,  the  corresponding  first-epoch  plate  area  was 
reprojected  onto  a  tangential  plane  centered  on  the  brick.  Then 
a  quadratic  solution  was  computed  to  transform  the  reprojected 
first-epoch  plate’s  (x,  y )  into  the  second-epoch  brick’s  (x,  y), 
forcing  the  chosen  reference  stars  to  have  zero  mean  proper 
motion.  The  quadratic  terms  in  the  solution  were  meant  to 
model  systematic  errors,  either  from  the  plate  (uncorrected 
Optical  Field  Angle  Distortion)  or  from  Tycho-2  proper-motion 
systematics10  unavoidably  propagated  backwards  in  time  into 
the  computed  first-epoch  (x,  y).  These  solutions  yield  measured 
relative  proper  motions  in  each  brick.  The  vast  majority  of  the 
solutions  only  needed  linear  terms  and  the  typical  standard  error 
of  the  solutions  varied  from  about  5  to  8  mas  yr_I,  which 
is  dominated  by  the  intrinsic  proper-motion  dispersion  of  the 
reference  stars. 

The  reference  stars  were  chosen  with  the  following  criteria: 
1  <  V  -  /  <  1.5, 0.25  <J-KS<  0.65  and  13  <  Ks  <  15, 
which  according  to  Nikolaev  &  Weinberg  (2000)  isolates  mostly 


10  Such  systematics  were  indeed  found  later  on  in  some  tests.  The  Tycho-2 
Catalog,  although  an  astrometric  catalog  based  on  late  epoch  space-based  data, 
has  early  epoch  positions  from  ground-based  data  that  are  known  to  suffer 
from  significant  magnitude  equation. 


G-M  dwarfs  in  the  Galaxy  disk,  located  between  0.4  to  1.6  kpc 
from  us,  with  an  estimated  mean  distance  of  «650  pc.  Their 
distribution  in  V  magnitude  ranges  mostly  from  V  =  15  to 
V  =  16.5.  In  general  there  are  between  200  and  500  reference 
stars  per  field,  depending  on  the  galactic  latitude  (our  fields 
extend  from  b  —50°  to  b  &  —20°).  The  intrinsic  proper- 
motion  dispersion  of  the  Galactic  disk  reference  stars  was  seen 
to  increase  with  \b\,  reflecting  the  changing  kinematics  of  the 
Galaxy  along  it. 

Despite  the  photometric  cuts  to  select  the  reference  stars, 
contamination  by  LMC  and  SMC  stars  could  not  be  avoided 
in  the  densest  parts  of  the  Clouds.  This  forced  us  for  the  time 
being  to  restrict  our  investigation  to  those  fields  in  which  we  trust 
the  relative  proper  motions,  as  being  measured  with  respect  to 
bona  fide  Galactic  foreground  stars.  These  areas  were  defined  as 
shown  in  Figure  3.  Although  a  substantial  number ofMagellanic 
Clouds  stars  were  lost  in  this  selection,  on  the  other  hand, 
confusion  due  to  image  crowding  at  these  locations  render  these 
fields  useless  anyway,  due  to  the  risk  of  misidentifications. 

From  the  initial  13880  bricks  available,  12180  are  in  the  non- 
contaminated  areas.  After  rejecting  bricks  from  the  discarded 
plates  (Section  3. 1),  10900  are  left  to  build  the  catalog  of  proper 
motions,  mostly  outside  the  Magellanic  Clouds.  In  order  to 
increase  the  number  of  LMC  and  SMC  stars  measured  at  the  end, 
contaminated  fields  that  overlap  with  this  non-contaminated 
catalog  were  later  on  directly  tied  into  it,  and  common  stars 
had  their  proper  motions  averaged. 

6.2.  Combining  the  Proper  Motions 

Once  relative  proper  motions  with  respect  to  the  Galactic 
foreground  stars  had  been  measured  in  the  non-contaminated 
area,  different  approaches  were  tried  to  combine  them  into  a 
single  global  well-defined  reference  frame.  Reference  stars  were 
chosen  hoping  that  their  mean  motion  along  the  sky  could  be 
described  by  a  smooth  function.  This  goal  was  indeed  attained  as 
confirmed  by  the  fact  that  the  measured  relative  proper  motion 
of  all  known  extragalactic  objects  were  very  precisely  fit  by  a 
quadratic  polynomial  in  (a,  5). 

At  this  point,  applying  this  polynomial  to  the  relative  proper 
motions  to  convert  them  to  absolute,  and  then  averaging  all 
measurements,  is  a  way  to  combine  all  the  information  available 
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Figure  4.  Systematics  when  comparing  a  single  frame  with  the  average  frame.  This  example  illustrates  the  case  of  a  PV  frame,  for  which  Am  =  MJ.singlc  frame  ■ 
ris.av cragc  frame  exhibits  visible  trends,  as  shown  by  the  median  line. 


per  star.  But  this  would  yield  a  catalog  with  a  rather  noisy  zero 
point  as  one  moves  along  the  sky,  mostly  due  to  the  real  intrinsic 
proper  motion  dispersion  of  the  reference  stars  and  the  fact  that 
two  overlapping  frames  may  not  have  the  same  reference  stars. 
Given  two  frames  with  about  50%  overlap  between  them,  means 
that  both  have  50%  of  the  reference  stars  in  common,  while 
the  other  50%  are  different.  Given  two  samples  with  N  data 
points  each,  both  with  the  same  dispersion  a  and  both  having 
50%  of  their  points  in  common,  it  can  be  easily  shown  that 

their  individual  mean  values  typically  differ  by  ■Jjjo.  In  other 
words,  we  can  expect  that  two  overlapping  frames  typically 
differ  in  their  relative  proper-motion  systems  by  «0.6  mas  yr~ 1 , 
for  N  —  220  and  a  =  6.5  mas  yr-1.  This  is  too  large  for  the 
level  of  precision  that  we  want  to  obtain. 

A  variant  of  the  so-called  block  adjustment  solutions 
(Eichhom  1960,  1988;  Jefferys  1979;  Stock  1981;  Taff  1988) 
was  considered  to  link  every  frame’s  relative  system  into  a  global 
one.  Each  frame,  containing  a  sample  of  the  whole  population 
of  reference  stars,  realizes  a  local  system  that  could  deviate 
from  the  global  system,  due  to  statistical,  measurement  and/ 
or  systematic  errors.  A  linear  function  per  frame  was  consid¬ 
ered  sufficient  to  describe  the  difference  between  the  frame  and 
the  global  system.  Therefore  three  parameters  per  frame,  per 
proper-motion  component  (jitt  cos  8  and  p.s  are  solved  sepa¬ 
rately)  need  to  be  determined.  The  coefficients  of  the  quadratic 
polynomial  that  globally  describes  the  mean  motion  of  the  ref¬ 
erence  stars  are  determined  as  well.  To  simultaneously  solve 
several  thousands  frames,  we  would  need  to  invert  an  approxi¬ 
mately  40,000  x  40,000  matrix  to  get  the  parameters  values  and 
their  errors.  For  this  project,  such  a  scheme  was  deemed  imprac¬ 
tical  at  this  time.  Nonetheless,  ideas  about  using  this  approach 
over  smaller  areas  first,  like  the  SPM  fields,  and  then  perform¬ 
ing  another  block  adjustment  solution  to  join  these  regions,  are 
being  adopted  for  future  work. 

A  more  practical  approach  was  considered  to  bring  each 
individual  frame’s  reference  system  closer  to  the  global  one. 
A  single  frame’s  reference  system  is  statistically  more  deviant 
from  the  global  system  than  is  the  system  defined  by  all  stars, 
since  the  proper-motion  reference  stars  are  only  a  subset  of  all 
stars.  Using  all  stars  in  common  between  overlapping  frames, 
we  can  adjust  each  frame’s  proper  motions  to  agree  with  the 
average  of  the  surrounding  fields.  This  adjustment  also  helps  to 
correct  residual  distortions,  as  they  are  statistically  smoothed 
out  in  the  average  frame.  To  avoid  frames  drifting  away  from 
the  global  system  as  they  are  being  aligned  to  one  another, 
all  reference  stars  in  the  field  are  explicitly  assigned  a  relative 
proper  motion  of  zero.  Once  the  adjustments  are  applied,  new 


averages  can  be  computed,  and  the  whole  process  is  iterated 
until  the  adjustments  converge  to  zero. 

Before  making  these  adjustments,  we  first  checked  to  see  if 
systematic  differences  existed  as  a  function  of  the  magnitude.  A 
non-negligible  number  of  frames  exhibited  systematic  trends 
with  magnitude  when  compared  to  the  average  frame  (see 
Figure  4).  In  general,  a  linear  function  of  the  coordinates  in 
the  field  would  take  care  of  the  geometrical  distortion,  but 
the  magnitude  equation  required  a  smoothed-localized  median, 
which  can  trace  the  general  trend  better  than  any  parameterized 
fit.  Therefore,  each  frame  is  first  corrected  for  its  differential 
magnitude  equation  with  respect  to  the  average  frame,  and  then 
we  proceed  to  correct  for  the  distortion,  following  the  iterative 
procedure  explained  above.  Each  frame  has  typically  about  1500 
to  2000  stars,  with  some  of  them  having  up  to  5000  stars,  to 
compute  the  adjustments. 

With  these  improved  relative  proper  motions,  the  quadratic 
polynomial  that  describes  the  mean  reflex  proper  motion  of 
the  extragalactic  objects  was  computed  and  used  to  transform 
the  relative  proper  motions  into  absolute  ones.  Polynomials 
for  fj.a  cos  5  and  p,s  as  functions  of  a  and  8  were  calculated 
separately,  using  proper-motion-error-weighted  least  squares. 
A  total  of  5351  external  galaxies  were  used  across  the  450 
sq-degree  area  and  the  formal  errors  of  the  polynomials,  com¬ 
puted  at  the  center  of  LMC  and  SMC,  based  on  the  full  co- 
variance  matrix,  amount  to  0.03  mas  yr-1  and  0.06  mas  yr-1, 
respectively. 

At  this  stage,  a  final  absolute  proper  motion  and  its  corre¬ 
sponding  error  are  obtained  per  star,  from  the  error-weighted 
average  of  all  the  individual  absolute  proper  motions  obtained 
for  it.  A  total  of  1 ,337,050  objects  in  the  non-contaminated  fields 
had  final  proper  motions,  including  a  good  number  of  LMC  and 
SMC  stars  located  in  the  outskirts  of  these  galaxies.  For  such 
reason,  it  is  named  the  outside  catalog. 

In  order  to  increase  the  number  of  LMC  and  SMC  stars 
with  measured  absolute  proper  motions,  contaminated  fields  that 
overlapped  with  the  outside  catalog  just  obtained,  were  directly 
tied  into  it  by  computing  a  linear  solution  to  correct  their  relative 
proper  motions  and  put  them  into  its  absolute  system.  A  total  of 
678  additional  fields  were  added  with  this  procedure,  eventually 
increasing  the  number  of  LMC  and  SMC  stars  by  about  30% 
and  50%,  respectively. 

6.3.  Zero-point  Global  Correction  of 
the  Absolute  Proper  Motions 

Since  we  are  now  theoretically  on  the  absolute  reference 
frame  defined  by  the  external  galaxies,  our  catalog  should  be 
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within  measurement  errors  in  the  same  reference  frame  system 
of  other  known  catalogs  of  absolute  proper  motion.  When  we 
checked  our  measures  (pa  cos  <5,  ps)  of  1356  Hipparcos  stars 
against  the  Hipparcos  Catalog,  we  found  a  significant  difference 
between  the  two  that  amounts  to 

Fa  cos  S, Hipparcos  ~  Fa  cos  S.This  work  =  0.49  i  0.07  mas  yr 

(1) 

and 

Fs,  Hipparcos  ~  Ms,  This  work  =  1.21  i  0.07  mas  yr  .  (2) 

The  source  of  this  systematic  difference  could  be  indeed  in 
any  (or  both)  of  the  two  catalogs.  Although  Hipparcos  is  the 
most  accurate  optical  astrometric  catalog  published  so  far,  it 
has  significant  correlations  between  the  astrometric  parameters 
(position,  proper  motion,  and  parallax)  of  different  stars,  when 
they  are  less  than  about  5°  apart  on  the  sky,  and  also  between 
the  astrometric  parameters  for  a  given  star,  due  to  the  special 
measurement  principle  of  Hipparcos  (Perryman  &  ESA  1 997). 
In  fact,  a  new  reduction  by  van  Leeuwen  (2007)  was  performed 
to  correct  some  systematic  correlations  in  the  data.  However, 
no  significant  difference  was  found  between  the  old  reference 
frame  (used  in  this  investigation)  and  the  new  one. 

On  the  other  hand,  our  plate  measurements  of  the  galaxies 
are  not  error-free.  Tests  were  run  to  compare  the  final  absolute 
proper  motions  if  only  galaxies  with  second-epoch  CCD  data,  or 
all  of  them,  were  used  to  compute  the  quadratic  polynomials  to 
transform  relative  proper  motions  into  absolute.  The  polynomial 
for  pu  cos  8  exhibited  noticeable  differences  around  the  SMC. 
This  was  not  completely  unexpected  given  that  at  that  location, 
four  SPM  fields  have  plate-only  data,  but  in  any  case  it  points  to 
the  fact  that  galaxies  in  the  second-epoch  plates  may  introduce 
problems.  Hence,  galaxies  in  the  first-epoch  plates  may  do  the 
same. 

As  determined  by  Girard  et  al.  (1998),  the  SPM  plate  material 
exhibits  magnitude  equation  for  galaxies  and  for  stars  that  differ 
in  functional  form.  However,  the  two  could  be  brought  into 
approximate  agreement  by  adding  an  offset  of  —0.7  to  the 
magnitudes  of  galaxy  images  before  calculating  the  magnitude- 
equation  correction  as  determined  from  stellar  images.  In  the 
present  study,  we  have  applied  the  same  —0.7  mag  offset  to 
galaxy  images  for  the  purpose  of  magnitude-equation  correction 
only. 

A  comparison  of  the  V  magnitudes  for  the  galaxies  in 
our  work,  showed  that  plate  photometry  returned  significantly 
fainter  magnitudes  than  in  CCD  calibrations.  Not  surprisingly, 
this  may  have  affected  the  magnitude  equation  correction  in  the 
plates,  since  it  was  the  calibrated  plate  photographic  magnitude 
that  was  used  for  such  purpose.  More  surprising  though,  was 
to  find  that  Equation  (1)  showed  a  linear  trend  versus  a  with  a 
slope  of  about  1 5  pas  yr- 1  per  degree  (a  bit  smaller  when  using 
only  galaxies  with  CCD  data).  The  LMC  and  SMC  center  of 
mass  positions  are  separated  by  about  20°,  for  which  the  trend 
above  indicates  a  zero  point  shift  of  0.3  mas  yr-1,  that  could 
in  principle  be  related  to  Hipparcos’  systematics.  On  the  other 
hand,  systematics  in  Equation  (2)  versus  5  were  also  seen,  that 
look  to  be  related  to  the  plates  location  and  layout. 

The  investigation  of  the  SPM  plates  magnitude  equation 
done  by  Girard  et  al.  (1998)  found  that  magnitude  equation 
terms11  varied  more  or  less  uniformly  from  first  to  second- 
epoch  plates,  so  we  can  expect  more  or  less  uniform  offsets 

1 1  A  polynomial  in  (X,  Y,  m)  is  used  to  describe  the  magnitude  equation 
correction. 


with  Hipparcos  proper  motions,  if  any  residual  uncorrected 
magnitude  equation  were  still  present  in  the  data.  Consequently, 
Equations  (1)  and  (2)  were  applied  to  all  the  absolute  proper 
motions,  putting  our  catalog  on  the  system  of  the  ICRS  via 
Hipparcos.  The  inaccuracy  of  its  zero  point  is  dominated  mostly 
by  Hipparcos ’  systematic  error  of  0.25  mas  yr-1,  since  the 
quadratic  polynomial,  as  defined  by  the  external  galaxies,  was 
in  general  very  accurate,  being  within  0.1  mas  yr-1  of  error  for 
most  of  the  450  deg2  area  studied. 

6.4.  Final  Catalog  of  Proper  Motions — Evaluation  of  Errors 

The  final  catalog  of  absolute  proper  motions  at  this  point 
has  1,448,438  objects,  with  the  following  data  listed:  a,  8,  V\ 
V  —  J,  J  —  Ks,  H  —  Ks  when  available;  absolute  proper  motions 
pa  cos  5,  ps  and  their  formal  errors  e/4|iC0S{,  s  in  mas  yr-1, 
number  of  data  points  used,  number  of  data  points  rejected 
(outliers  were  rejected  based  on  their  normalized  errors),  a 
flag  to  indicate  Hipparcos,  Tycho-2,  2MASS  extended  sources, 
confirmed  LEDA  Galaxies  and  QSOs,  a  flag  to  indicate  if  the 
object  is  or  is  not  a  reference  star,  and  the  2MASS  ID  when 
available.  The  overall  distribution  of  the  stars  in  the  catalog  can 
be  seen  in  Figure  5. 

In  the  final  catalog  of  proper  motions,  stars  with  V  <  12 
have  formal  proper  motion  errors  of  about  0.5  mas  yr-1,  and 
well-measured  stars  with  12  <  V  <  15.5,  have  values  that 
range  from  0.5  to  1.3  mas  yr-1 .  These  only  reflect  measurement 
errors,  as  they  are  based  on  the  positional  errors  and  the  epoch 
difference.  When  proper  motions  for  the  same  star  measured  in 
different  bricks  are  averaged  together,  statistical  and  systematic 
deviations  between  the  bricks  must  be  added  to  get  the  true 
error.  The  iterative  method  applied  in  Section  6.2  was  designed 
to  reduce  those  deviations,  but  cannot  make  them  zero. 

A  better  way  to  determine  the  real  proper  motion  uncertain¬ 
ties,  is  to  compute  the  standard  deviation  of  the  final  error- 
weighted  average  proper  motions.  Figure  6  shows  these  scatter- 
based  proper  motion  errors,  for  a  random  sample  of  about  5% 
the  size  of  the  whole  catalog.  These  errors  can  still  be  an  un¬ 
derestimate  of  the  real  proper  motion  uncertainty,  because  in 
the  error-weighted  average  the  data  points  are  not  independent, 
two  overlapping  bricks  reduced  into  the  same  plate  can  pro¬ 
duce  (positively)  correlated  proper  motions.  In  this  case,  the 
true  uncertainty  is  larger  than  the  measured  one. 

The  most  reliable  assessment  of  the  proper  motion  errors  is 
obtained  from  a  comparison  with  external  catalogs.  The  scatter 
observed  in  the  differences  between  our  proper  motions  and 
those  from  Hipparcos,  is  the  combined  result  of  both  catalogs’ 
proper  motion  errors.  Given  that  Hipparcos  errors  are  about 
1  mas  yr-1,  the  measured  dispersion  indicates  that  our  real 
proper  motion  uncertainties  are  about  2.3  mas  yr-1,  for  stars 
brighter  than  V  =  10.  This  coincides  well  with  the  scatter- 
based  errors  in  Figure  6,  at  the  bright  end. 

The  scatter  in  the  proper  motion  of  the  LMC  and  SMC 
samples  can  be  used  to  estimate  the  proper  motion  uncertainties, 
at  the  mean  magnitude  of  the  Clouds,  since  their  intrinsic  internal 
velocity  dispersion12  makes  a  very  minor  contribution  to  the 
observed  scatter.  Results  indicate  that  our  real  proper  motion 
uncertainties  are  about  3.8  mas  yr-1  for  stars  around  V  = 
16.4,  entirely  consistent  with  the  scatter-based  uncertainties  in 
Figure  6. 


12  30  km  s  1  (0.13  mas  yr  1  at  50  kpc  distance)  or  less  in  each  (Westerlund 
1997). 
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Figure  5.  Sky  coverage  of  the  final  catalog  of  proper  motions.  Every  tenth  entry  is  plotted  in  light  gray.  The  “contaminated”  areas  are  as  in  Figure  3.  The  lower  density 
area  north  of  the  SMC  corresponds  to  an  area  where  several  poor  quality  plates  were  discarded. 


V  (mag) 

Figure  6.  Uncertainty  in  in  vs.  magnitude.  The  error  estimate  is  based  on 
the  scatter  of  the  individual  data  points  that  were  error-weighted-averaged  into 
the  final  proper  motion.  The  black  line  is  a  moving  median,  computed  over  a 
0.5  mag  interval.  The  gray  dashed  line  is  the  moving  median  value  for  u<t  cos  5. 
This  plot  is  based  on  a  random  sample  of  5%  the  size  of  the  whole  catalog. 

At  the  faint  end,  the  dispersion  in  the  proper  motion  of 
external  galaxies  varies  substantially,  depending  on  whether 
they  have  CCD  or  plate  second-epoch  data.  In  the  first  case, 
the  dispersion  is  1 1  mas  yr_I  at  a  mean  magnitude  of  V  =  17.3. 
In  the  second  case,  the  dispersion  is  21  mas  yr-1  at  a  mean 
magnitude  of  V  =  18.5.  Since  the  plate  photometry  for  the 
galaxies  produced  systematically  fainter  magnitudes  than  the 
CCD  photometry,  the  difference  in  these  two  sets  reflects 
also  the  difference  in  precision  between  (better)  CCD  and 
(worse)  plate  measurements,  and  not  just  the  increase  of  errors 
with  magnitude.  At  these  magnitudes,  our  scatter-based  proper 
motion  uncertainties  are  about  50%  below  these  values. 

In  general,  the  scatter-based  proper  motion  errors  are  a  rather 
good  indicator  of  the  real  uncertainties  in  the  proper  motions 
for  stars  brighter  than  about  V  =  16.5.  More  importantly, 
all  the  above  external  estimates:  2.5,  3.8  and  11  mas  yr-1  at 
V  ~  10, 16.4  and  17.3,  respectively,  are  smaller  in  size  than 
what  was  achieved  in  SPM3,  clearly  showing  the  increased 
precision  due  to  having  second-epoch  CCD  data. 


7.  PROPER  MOTION  OF  THE  MAGELLANIC  CLOUDS 
7.1.  Selection  ofLMC  and  SMC  Dominated  Samples 

A  photometric  selection  was  made  to  choose  bona  fide  red 
giant  LMC  and  SMC  stars,  based  on  the  analysis  of  the  2MASS 
LMC  infrared  color  magnitude  diagram  (CMD)  of  Nikolaev 
&  Weinberg  (2000)  (their  sample  “J”).  The  photometric  cuts 
applied,  as  seen  in  Figure  7,  are 

LMC:  1.1  s?  J  -  Ks  <  1.3  and  9.5  sj  Ks  St  12, 
SMC:  1.0  <  J  -  Ks  <  1.2  and  10  <  Ks  <  12.5. 

Only  stars  with  CCD  second-epoch  data  were  selected,  as  stars 
whose  proper  motions  were  based  on  plate  data  only  showed 
a  significantly  higher  dispersion  and  some  visible  systematics. 
We  were  also  forced  to  discard  an  area  of  high  stellar  density, 
at  71°  <  a  <  76°  and  -68°  ^  S  $  —71°,  with  CCD  data  close 
to  the  LMC  center,  which  consistently  showed  deviant  results 
probably  caused  by  misidentifications. 

We  selected  3822  LMC  and  964  SMC  stars,  as  seen  in 
Figure  8,  to  measure  the  mean  absolute  proper  motion  of  the 
Clouds.  Bluer  sequences  of  the  CMDs  in  Figure  7  containing 
LMC/SMC  stars  have  significant  contribution  from  Milky  Way 
stars  and  therefore  were  not  considered.  The  redder  sequences 
of  LMC/SMC  asymptotic  giant  branch  stars  have  very  faint 
V  «  17.5  mag,  consequently  the  proper  motion  errors  are  too 
large  to  be  useful.  Also,  given  the  magnitude-related  problems  in 
the  plates,  it  is  desirable  to  have  the  smallest  possible  difference 
in  brightness  between  the  reference  stars  and  the  Clouds’ 
stars,  and  indeed,  the  chosen  samples  overlap  sufficiently  in 
magnitude  (see  Figure  9). 

7.2.  Absolute  Proper  Motion  of  the  LMC  and  SMC 

Probability  plots  (Hamaker  1978)  of  the  chosen  samples 
yielded  the  mean  and  dispersion  values  for  the  LMC  and  SMC 
proper  motion  listed  in  Table  2.  The  errors  quoted  include: 
the  formal  error  of  the  mean  value  {cr/^/Nsws),  the  error 
of  the  quadratic  polynomial  at  the  LMC  and  SMC  centers, 
transformation  to  Hipparcos  errors  (Equations  (1)  and  (2)),  and 
the  estimated  Hipparcos  systematic  error  (0.25  mas  yr-1).  As 
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Figure  7.  Infrared  CMDs  in  the  vicinity  of  the  LMC  (left)  and  SMC  (right)  based  on  2MASS  photometry  of  the  stare  in  our  proper  motion  catalog.  The  two  bluest 
vertical  sequences  visible  in  both  CMDs  are  dominated  by  foreground  Galactic  stars,  while  the  bright  red  sequences  are  dominated  by  LMC/SMC  stars.  Black  points 
on  each  panel  show  the  red  giant  star  samples  selected  to  measure  the  mean  motion  of  each  Cloud. 


Figure  8.  Distribution  on  the  sky  of  LMC  and  SMC  sample  stars  selected  from  our  catalog  of  proper  motions. 


v 


Figure  9.  Histogram  of  Vmagnitude  for  the  reference  stars,  LMC  and  SMC  stars. 
All  three  data  sets  overlap  substantially  in  V,  which  helps  to  reduce  systematics 
in  the  proper  motions  related  to  brightness. 

explained  before,  the  error  budget  is  dominated  by  Hipparcos 
systematics. 

Table  3  and  Figures  10  and  i  1  summarize  how  our  results 
compare  with  recent  measurements  of  the  proper  motion  of 
the  Magellanic  Clouds.  The  error  bars  hinder  a  more  precise 


t  Reference  stars 
f!  SMC  stars 
ft  LMC  stars 


conclusion  about  the  individual  tangential  velocities  of  the 
Clouds  based  on  our  data.  Nonetheless,  the  methodology  used 
to  measure  the  stellar  proper  motions  in  our  catalog  permits  us 
to  make  a  rather  precise  measurement  of  the  proper  motion  of 
the  SMC  with  respect  to  the  LMC,  as  explained  in  more  detail 
in  Section  7.3. 


7.2.1.  Center-of-mass  Proper  Motion 

The  large  extent  of  the  Clouds  and  their  non-negligible  depth 
means  that  all  previous  investigations,  which  measured  proper 
motions  on  scattered  small  fields,  had  to  convert  their  measured 
values  into  a  center-of-mass  proper  motion.  That  is  because  a 
given  space  velocity  at  a  fixed  distance  projects  differently  on 
radial  velocity  and  proper  motion  at  different  locations  in  the 
sky,  following  the  same  principle  of  the  moving  cluster  method 
and  the  solar  motion.13  Besides,  proper  motion  obviously  scales 
with  distance. 

13  One  of  the  first  attempts  to  estimate  the  proper  motion  of  the  Magellanic 
Clouds  was  done  by  measuring  gradients  in  the  radial  velocity  along  them 
(Feitzinger  et  al.  1977;  Mealheringham  et  al.  1988). 
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Table  2 

Absolute  Proper  Motion  of  the  Magellanic  Clouds 


Target 

■Wstare 

COS  8 

(mas  yr-1) 

Mi 

(mas  yr-1) 

Zfla  COS  8 

(mas  yr-1) 

(masyr  ') 

Gfla  COS  8 

(mas  yr-1) 

(mas  yr-1) 

LMC 

3822 

1.89 

0.39 

0.27 

0.27 

3.76 

3.59 

SMC 

964 

0.98 

-1.01 

0.30 

0.29 

4.12 

3.82 

Notes.  The  e  values  are  the  errors  of  the  mean  p  values,  which  include:  the  formal  errors  (_<? / Nsi3rj),  the  error  of  the 
quadratic  polynomial  at  the  LMC  and  SCM  centers,  transformation  to  Hipparcos  errors  (Equations  (1)  and  (2)),  and  the 
Hipparcos  systematic  error  (0.25  mas  yr-1 ).  The  a  values  represent  the  scatter  of  the  data  around  the  mean  proper  motion 
values. 


Figure  10.  Proper  motion  of  the  Large  Magellanic  Cloud.  Labels  are  as  follows: 
J94,  Jones  etal.(  1994);  K97,  Kroupa  &  Bastian  (1997c);  K06,  Kallivayalil  et  al. 
(2006b);  P06,  Pedreros  et  al.  (2006a);  P08,  Piatek  et  al.  (2008);  C09-1,  Costa 
et  al.  (2009)-(1)  and  C09-2,  Costa  et  al.  (2009)-(2).  All  these  results  are  plotted 
in  gray  while  ours  is  plotted  in  black. 

Additionally,  internal  rotation  must  also  be  taken  into  account 
for  the  LMC.  In  the  case  of  the  SMC,  such  correction  is  deemed 
unnecessary  because  its  stellar  component  is  mostly  supported 
by  velocity  dispersion  (Harris  &  Zaritsky  2006).  The  LMC’s 
rotation  curve  is  obtained  from  radial  velocities  of  Carbon 
stars  (van  der  Marel  et  al.  2002;  Olsen  &  Massey  2007)  and 
yield  widely  accepted  values  of  Vrot.LMc  =  50-60  km  s-1. 
Nonetheless,  Piatek  et  al.  (2008)  estimates  its  own  Vro!j,MC  = 
120  km  s_1,  based  on  the  gradient  of  their  measured  proper 
motions  along  the  radius  in  the  LMC  disk,  and  use  such  value 
for  the  rotation  correction. 

Given  such  discrepant  values,  Costa  et  al.  (2009)  actually 
obtains  two  final  results  for  the  LMC,  Costa  et  al.  (2009)- 
(1)  refers  to  their  final  proper  motion  when  using  V^mc  = 
50  km  s-1,  while  Costa  et  al.  (2009)-(2)  does  so  for  Vrot,LMC  = 
120  km  s-1 .  For  the  latter,  it  must  be  noted  that  cos  S  deviates 
noticeably  from  the  other  determinations  in  Table  3. 


_ i _ i _ i _ i _ . _ i _ i _ i_ 
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Proper  Motion  RA*  Cos(DEC)  {mas/yrj 

Figure  II.  Proper  motion  of  the  Small  Magellanic  Cloud.  Labels  are  as  follows: 
K97,  Kroupa  &  Bastian  (1997c),  K06,  Kallivayalil  et  al.  (2006a),  P08,  Piatek 
et  al.  (2008),  and  C09,  Costa  et  al.  (2009).  All  these  results  are  plotted  in  gray 
while  ours  is  plotted  in  black. 

Table  3 

Recent  Determinations  of  the  Proper  Motion  of  the  Magellanic  Clouds 


Author  pa  cos  8  ps  JVaars 

_ (mas  yr~*) _ (mas  yr-1) _ 


Large  Magellanic  Cloud 


Jones  etal.  (1994) 

1.37 

0.18 

0.28 

£ 

0.27 

251 

Kroupa  &  Bastian  (1997c) 

1.94 

± 

0.29 

-0.14 

± 

0.36 

33 

Pedreros  et  al.  (2006a) 

1.80 

xfc 

0.10 

0.90 

£ 

0.10 

108 

Kallivayalil  et  al.  (2006b) 

2.03 

£ 

0.08 

0.44 

£ 

0.05 

810 

Piatek  et  al.  (2008) 

1.95 

± 

0.04 

0.44 

± 

0.04 

889 

Costa  etal.  (2009)  (1) 

1.82 

± 

0.13 

0.39 

± 

0.15 

41 

Costa  etal.  (2009)  (2) 

1.61 

± 

0.13 

0.60 

£ 

0.15 

41 

This  work 

1.89 

£ 

0.27 

0.39 

± 

0.27 

3822 

Small  Magellanic 

Cloud 

Kroupa  &  Bastian  (1997c) 

1.23 

£ 

0.84 

-1.21 

± 

0.75 

9 

Kallivayalil  et  al.  (2006a) 

1.16 

£ 

0.18 

-1.17 

£ 

0.18 

177 

Piatek  et  al.  (2008) 

0.75 

£ 

0.06 

-1.25 

± 

0.06 

215 

Costa  et  al.  (2009) 

1.03 

± 

0.29 

-1.09 

£ 

0.18 

44 

This  work 

0.98 

± 

0.30 

-1.10 

± 

0.-29 

964 

Notes.  Results  (1)  and  (2)  from  Costa  et  al.  (2009)  are  obtained  assuming  two 
different  rotational  velocities  for  the  LMC,  as  explained  in  detail  in  Section  7.3. 

Altogether,  the  typical  correction  for  perspective  effect  for 
the  LMC  from  both  methods  is  about  ±  0.2  mas  yr-1  and  may 
rise  to  about  ±  0.5  mas  yr-1  for  locations  farther  than  5°  from 
the  LMC  center,  running  more  or  less  in  opposite  directions  at 
opposite  locations  on  the  Cloud.  Rotation  effects  for  the  LMC 
are  usually  less  than  0.1  mas  yr-1.  Perspective  effects  for  SMC 
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Figure  12.  Comparison  of  previous  and  new  determinations  of  the  Proper  motion 
of  the  SMC.  Previous  determinations  from  Figure  1 1  are  plotted  with  filled 
translucent  symbols:  K97  in  yellow,  K06  in  blue,  P08  in  green,  and  C09  in  red. 
Our  new  determinations  based  on  the  same  studies’  LMC  proper  motion  plus 
our  relative  SMC-LMC  proper  motion,  are  plotted  with  open  symbols  of  the 
same  colors  (C09  has  two  determinations  of  the  LMC  proper  motion).  Previous 
results  are  consistent  with  ours  when  the  filled  and  the  open  symbol  of  the  same 
color  overlap.  Our  SMC  proper  motion  in  plotted  with  the  black  symbol. 


are  much  smaller.  For  previous  studies,  these  corrections  are 
necessary  and  in  some  cases  yield  quite  different  values  for 
the  same  fields,  but  in  our  work,  given  the  spatial  extent  and 
symmetry  of  the  data,  the  net  effect  on  the  mean  motion  of  the 
Clouds  is  very  close  to  zero,  and  no  correction  is  done. 

7.3.  Relative  Proper  Motion  of  the  SMC 
with  Respect  to  the  LMC 

As  explained  in  Section  6.2,  we  measured  the  mean  motion 
of  our  reference  stars  precisely  all  over  our  FOV,  in  particular  at 
the  location  of  LMC  and  SMC  within  0.03  and  0.06  mas  yr-1, 
respectively.  Combined  with  the  relative  proper  motion  of  LMC 
and  SMC  stars  with  respect  to  these  reference  stars,  we  can 
indeed  measure  the  proper  motion  of  the  SMC  with  respect  to 
that  of  the  LMC,  with  a  higher  precision,  limited  by  the  error 
just  quoted  plus  the  formal  error  of  the  mean  coming  from  the 
number  of  stars  and  their  measured  scatter. 

From  Table  2  and  taking  the  errors  quoted  above  into  account, 
it  is  straightforward  to  obtain  the  relative  proper  motion  of  the 
SMC  with  respect  to  the  LMC: 

&pa cos s (SMC-LMC)  =  -0.91  ±  0.16  mas  yr-1  (3) 

Ams(SMC-LMC)  =  -  1.49  ±  0.15  mas  yr-1.  (4) 

These  values  cannot  be  transformed  directly  into  a  measurement 
of  the  relative  velocity  between  the  Clouds.  Being  at  different 
locations  in  the  sky  means  that  the  planes  of  their  tangential 
velocities  are  different  as  well,  and  the  necessary  rotation  and 
projections  to  measure  the  SMC  velocity  on  the  LMC  reference 
frame  does  not  allow  us  to  obtain  the  relative  space  velocity 
as  merely  a  function  of  the  relative  proper  motion  between  the 
Clouds. 


Table  4 

New  Determinations  of  the  Proper  Motion  of  the  SMC 


Using  the  LMC 

Proper  Motion  of 

cos  & 

(mas  yr-1) 

US 

(mas  yr' 

-') 

Jones  et  al.  (1994) 

0.46 

± 

0.31 

-1.21 

± 

0.31 

Kroupa  &  Bastian  ( 1 997c) 

i.03 

± 

0.39 

-1.63 

0.39 

Pedreros  et  al.  (2006a) 

0.89 

± 

0.18 

-0.59 

± 

0.19 

Kallivayalil  et  al.  (2006b) 

1.12 

± 

0.16 

-1.05 

0.18 

Piatek  et  al.  (2008) 

1.04 

± 

0.16 

-1.05 

± 

0.16 

Costa  et  al.  (2009)  (1) 

0.90 

0.21 

-1.10 

± 

0.21 

Costa  et  al.  (2009)  (2) 

0.70 

± 

0.21 

-0.60 

± 

0.21 

Notes.  Results  (1)  and  (2)  from  Costa  et  al.  (2009)  are  obtained  assuming  two 
different  rotational  velocities  for  the  LMC,  as  explained  in  detail  in  Section  7.3. 
Compare  with  the  direct  determinations  of  the  proper  motion  of  the  SMC  in 
Table  3.  Quoted  errors  include  the  contribution  from  our  measured  relative 
proper  motion. 

But,  we  can  use  these  values  to  obtain  new  independent  mea¬ 
surements  of  the  SMC’s  absolute  proper  motion,  based  on  ex¬ 
isting  measurements  of  the  LMC’s  absolute  proper  motion  plus 
our  precise  relative  proper  motion  from  above.  Moreover,  since 
all  authors  that  directly  measured  the  proper  motion  of  the  SMC 
had  previously  measured  the  LMC’s  proper  motion  as  well,  we 
can  verify  if  their  original  SMC  results  are  consistent  with  our 
relative  measure.  Figure  12  shows  the  absolute  proper-motion 
determination  for  the  SMC  from  Table  4  (made  by  combining 
our  relative  LMC-SMC  motion  with  absolute  LMC  motions 
from  the  literature),  compared  to  the  direct  determinations  of 
the  absolute  proper  motion  of  the  SMC  from  Table  3. 

Except  for  Piatek  et  al.  (2008)  and  Costa  et  al.  (2009)-(2), 
all  published  measurements  of  the  SMC  proper  motions  are 
consistent  with  our  new  ones.  These  two  works  are  the  only 
ones  that  use  VrocLMC  =  120  km  s-1.  It  is  worth  noting  that 
field  LI  1  in  the  LMC  from  Piatek  et  al.  (2008)  is  also  in  Costa 
et  al.  (2009),  but  their  measured  proper  motions  are  significantly 
different,  beyond  the  quoted  errors.  All  this  makes  us  suspect 
that  the  rather  small  quoted  proper  motion  errors  in  Piatek  et  al. 
(2008)  underestimate  their  real  uncertainties. 

In  summary,  although  our  proper  motions  for  the  LMC  and 
the  SMC  separately  are  in  agreement — within  error  bars — with 
Kallivayalil  et  al.  (2006b,  2006a),  Piatek  et  al.  (2008),  and  Costa 
et  al.  (2009),  our  relative  proper  motion  of  SMC  with  respect  to 
LMC  is  consistent  only  with  Kallivayalil  et  al.  (2006a,  2006b) 
and  Costa  et  al.  (2009)-(l). 

7.4.  The  Space  Motion  of  the  Clouds 

The  individual  as  well  as  relative  space  velocities  of  the 
Clouds,  as  derived  from  our  LMC  and  SMC  proper  motions, 
are  given  in  Table  5.  The  escape  velocity  at  the  distance  of 
the  LMC  is  estimated  to  be  300-350  km  s“',  depending  on 
the  Galactic  potential  model  used  (either  a  simple  isother¬ 
mal  sphere  or  a  more  elaborate  “cosmologically  inspired” 
Navarro-Frank- White  (NFW)  dark  matter  profile  (Gardiner  & 
Noguchi  1996;  Klypin  et  al.  2002;  Xue  et  al.  2008).  Taken  at 
face  value,  the  galactocentric  velocities  in  Table  5  indicate  that 
the  LMC  is  traveling  at  a  speed  that  is  very  close  to  the  escape 
velocity,  while  the  SMC  is  still  below  the  escape  velocity  of  the 
Galaxy.  Unfortunately,  the  uncertainties  hinder  a  more  definitive 
conclusion  regarding  their  binding  status. 

From  our  proper  motion  measures,  we  determine  that  the 
SMC  is  moving  at  89  ±  54  km  s-1  with  respect  to  LMC.  Our 
error  bars  do  not  allow  us  to  determine  whether  or  not  the  Clouds 
are  bound  to  each  other,  given  that  the  escape  velocity  from  the 
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Tables 

Space  Velocity  of  the  Magellanic  Clouds 


Parameter 

Unit 

LMC 

SMC 

Measured  or  Adopted  Quantities 

(«.«) 

deg 

(81.90,  -69.87) 

(13.20,  -72.50) 

(Lb) 

deg 

(280.526,  -32.527) 

(303.788,  -44.628) 

Distance 

kpc 

50.1  ±2.3 

62.8  ±  2.6 

Radial  velocity 

km  s~l 

262.1  ±3.4 

146.0  ±  0.6 

OiaCOS 8,  m) 

masyr-1 

(1.89, 0.39) 

(0.98,-1.10) 

±  (0.27,  0.27) 

±  (0.30, 0.29) 

( U .  V,  HOo.lsr 

km  s"1 

(10.00,  5.25, 7.17) 

(U,  V.  W)lsr.8c 

km  s-1 

(0, 220, 0) 

Derived  Quantities 

(X,  Y,  Z) 

kpc 

(0.3,  -41.5,  -26.9) 

(-16.2,  -37.6,  -44.1) 

Rgc 

kpc 

41.5 

40.9 

(U,  V,  W)hc 

km  s~‘ 

(71.8,  -468.1,226.1) 

(51.5,  -433.4, 139.2) 

±(63.6,37.1,56.6) 

±(82.1,70.3,62.2) 

(U,  V,  W)sc 

km  s_l 

(71.8,  -248.1,226.1) 

(51.5,  -213.4, 139.2) 

±(63.6,37.1,56.6) 

±(82.1,70.3,62.2) 

lice/,  V,  W)sc|l 

km  s-1 

343.3  ±  47.8 

259.9  ±  68.6 

ll(f/.  v.  W)smc-lmc!I 

km  s_1 

89.2  ±  53.9 

Notes,  (a,  S )  from  van  der  Marel  et  al.  (2002)  (LMC)  and  Stanimirovic  et  al.  (2004)  (SMC).  (I,b)  are 
transformed  from  (cr,  5).  Distance  from  Freedman  et  al.  (2001)  (LMC)  and  Cioni  et  al.  (2000)  (SMC). 
Radial  velocity  from  van  der  Marel  et  al.  (2002)  (LMC)  and  Harris  &  Zaritsky  (2006)  (SMC).  Solar 
motion  ((/,  V,  W)q,lsr,  and  galactocentric  velocity  of  the  Local  Standard  of  Rest  ( U ,  V,  lV)LSR,gC.  from 
Dehnen  &  Binney  (1998).  R0  =  8  kpc  was  adopted.  Rgc  is  the  Galactocentric  distance.  ( U ,  V,  IV);, c  and 
( U .  V,  W)gc  are  the  heliocentric  and  galactocentric  space  velocities  of  the  Clouds.  ||((/,  V,  W)lmc,smcII 
is  the  relative  velocity  of  SMC  with  respect  to  LMC. 


Table  6 

Relative  Velocity  between  LMC  and  SMC 


Reference 

ll(t/,  V,  W)SMC-lmcII 
km  s-1 

Originally  quoted  by  reference 

Kallivayalil  et  al.  (2006a) 

105  ±  42 

Piatek  et  al.  (2008) 

142  ±  19 

Costa  et  al.  (2009)-(l) 

84  ±  50 

Costa  et  al.  (2009)-(2) 

62  ±  63 

New  value  using  our  fiSMC-LMC 

Kallivayalil  et  al.  (2006a) 

123  ±  51 

Piatek  et  al.  (200S) 

108  ±  51 

Costa  et  al.  (2009)-(l) 

83  ±  52 

Costa  et  al.  (2009)-(2) 

152  ±  63 

This  work 

This  work 

89  ±  54 

Notes.  Results  (1)  and  (2)  from  (Costa  et  al.  2009)  are 
obtained  assuming  two  different  rotational  velocities  for 
LMC,  as  explained  in  detail  in  Section  7.3. 


LMC  at  the  SMC  location  is  about  90  km  s_1  (assuming  a  simple 
point-mass  geometry,  a  mass  for  the  LMC  of  2  x  1010  MQ  and 
23  kpc  for  the  distance  between  the  Clouds).  Yet,  we  can  use  our 
new  SMC  proper  motions  to  obtain  additional  determinations 
of  the  relative  velocity  between  the  Clouds,  using  the  more 
precise  LMC  proper  motions  available  in  the  literature,  and  the 
other  needed  input  parameters  listed  in  Table  5.  Then  we  obtain 
more  estimates  of  ||((7,  V,  W)SMC_LMcll>  which  are  listed  in 
Table  6.  For  comparison,  this  table  also  lists  the  relative  velocity 
originally  quoted  by  the  references  used. 

It  is  important  to  be  aware  of  how  sensitive 
1 1  (U,  V,  W)SMC_LMC  1 1  is  for  different  values  of  the  LMC  proper 


motion  (see  Figure  13).  It  is  clear  that  even  the  smallest  proper 
motion  errors  quoted  translate  into  a  substantial  relative  ve¬ 
locity  error  between  the  Clouds.  On  top  of  that,  the  mea¬ 
sured  values  are  close  enough  to  the  escape  velocity  of  the 
SMC  with  respect  to  the  LMC,  such  that  any  conclusion  re¬ 
garding  the  binarity  of  the  Clouds  is  still  far  from  decided. 
The  same  situation  applies  to  the  individual  space  veloci¬ 
ties  of  the  Clouds,  the  binding  status  of  the  Clouds  to  the 
Milky  Way  is  also  extremely  sensitive  to  the  individual  proper 
motions. 

Thus,  although  our  individual  absolute  proper  motions  for 
the  Clouds  cannot  tell  us  much  about  their  orbits,  our  measured 
relative  proper  motion  between  LMC  and  SMC  allowed  us  to 
identify  the  best  self-consistent  measurements  of  the  proper 
motion  of  both  clouds,  from  those  authors  that  had  measured 
both  galaxies  with  claimed  very  good  accuracies.  Our  measured 
space  velocity  has  also  helped  to  check  which  results  were  more 
consistent  with  that  value,  since  it  is  quite  sensitive  to  the  proper 
motions  used. 

8.  IMPLICATIONS  OF  OUR  RESULTS  TO  THE 

DYNAMICS  OF  THE  MAGELLANIC  SYSTEM 

The  most  prevalent  scenario  for  the  orbit  of  the  Magellanic 
Clouds,  before  the  HST  proper  motion  results,  favored  them 
as  a  binary  system  in  a  bound  orbit  around  the  Milky  Way. 
Dynamical  simulations  were  performed  under  the  presumption 
that  the  Clouds  were  bound  to  each  other,  due  to  the  existence 
of  the  common  H  i  envelope  surrounding  them  (Hindman  et  al. 
1 963).  Crude  timing  estimates  suggested  that  for  its  creation  and 
survival,  a  long  time  of  shared  orbits  was  needed.  The  extent  of 
the  Magellanic  Stream  implied  that  the  Clouds  had  undergone 
multiple  orbits  gravitationally  bound  to  the  Milky  Way  and  to 
each  other. 
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Figure  13.  Relative  velocity  between  the  Clouds  as  a  function  of  the  absolute 
proper  motion  of  the  LMC.  For  any  given  value  ()ia  cosS,  us)lmc  in  the 
diagram.  cosS.  us)smc  is  obtained  from  Equations  (3)  and  (4)  and  then 
ll(J/,  V.  HOsmc— lmcII  is  computed  (other  needed  parameters  from  Table  5) 
and  plotted  with  a  color  (gray  in  the  printed  journal)  indicative  of  its  value. 
LMC  proper  motions  and  errors  from  K06,  P08,  C09-1,  and  C09-2  are  plotted 
in  black  with  their  error  ellipses,  and  from  this  work  with  error  bars  only.  The 
black  parabolic  curve  indicates  the  values  of  LMC  proper  motion  for  which 
||(t/,  V,  WOsmc-LMcII  =  90  km  s-1.  Values  of  LMC  proper  motion  that  fall 
below  this  curve  roughly  correspond  to  the  Clouds  forming  a  bound  binary 
system. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

The  orbits  of  the  Clouds  are  naturally  a  key  ingredient  in 
such  dynamical  modeling.  Once  an  orbit  is  chosen,  then  a  full 
simulation  of  the  Clouds  themselves  is  run.  The  main  goal  has 
been  to  reproduce  as  much  as  possible  the  position  and  radial 
velocity  of  the  Magellanic  Stream,  though  significant  validation 
is  obtained  if  the  Leading  Arm,  the  Magellanic  Bridge,  and  the 
distorted  structure  of  the  Clouds,  are  also  replicated.  Another 
“condition”  imposed  is  that  the  Clouds  share  the  same  orbital 
plane  with  the  Magellanic  Stream,  since  the  latter  runs  more  or 
less  on  a  great  circle  on  the  sky. 

Depending  on  the  particular  questions  under  study,  Cloud 
model  investigators  use  different  approaches.  Either  one  or 
both  Clouds  are  considered,  using  collisionless  (stars)  and/or 
collisional  (gas)  particles.  The  Cloud(s)  is(are)  made  of  one 
or  many  particles,  which  can  be  massless  test  particles  under  a 
given  potential,  or  a  conglomerate  of  self-gravitating  particles  or 
even  “sticky”  self-gravitating  particles  (to  model  hydrodynamic 
processes). 

Overall,  two  competing  scenarios  have  been  systematically 
studied  and  have  had  some  degree  of  success  in  the  modeling 
of  the  mentioned  structures:  the  tidal  model  and  the  ram 
pressure  model.  In  the  first  scenario,  the  Galaxy  extracted  a 
tidal  plume  from  the  LMC  and/or  the  SMC,  which  gave  origin  to 
the  Magellanic  Stream,  in  a  previous  close  encounter  of  the  three 
bodies  about  1 .5  Gyr  ago.  A  most  recent  encounter  200  Myr  ago 
created  the  Magellanic  Bridge.  The  second  scenario  proposes 
that  the  Stream  and  the  Leading  Arm  consist  of  material  that  has 
been  ram  pressure  stripped  from  the  LMC  (and  SMC),  during  its 
last  passage  through  the  extended  ionized  Halo  of  the  Galaxy, 
about  500  Myr  ago. 

A  complementary  scenario  to  extract  substantial  amounts  of 
gas  and  no  stars  from  the  Clouds,  the  blowout  model,  has  been 
proposed  by  Nidever  et  al.  (2008).  The  main  hypothesis  is  that 
star  formation  outflow  in  the  leading  edge  of  the  LMC  has  been 


blowing  out  or  puffing  up  the  gas  over  the  past  2  Gyr,  making  it 
easier  for  ram  pressure  and  tidal  forces  to  strip  it  off. 

The  first  dynamical  simulation  of  the  Magellanic  System 
(Murai  &  Fujimoto  1980)  was  done  knowing  only  the  radial 
velocity  of  the  Clouds,  and  very  little  about  the  mass  of 
the  Milky  Way  at  the  Clouds’  distances.  Some  dynamically 
permissible  parameters14  were  assumed  in  order  to  have  all 
the  input  parameters  necessary  to  compute  an  orbit.  Their 
models  required  the  inclusion  of  a  massive  halo  (an  idea  that 
was  just  starting  to  being  accepted  then)  in  order  to  reproduce 
the  highly  negative  radial  velocities  observed  in  the  tip  of  the 
Magellanic  Stream.  Conversely,  assuming  that  the  LMC  was 
in  a  bound  orbit  provided  estimates  consistent  with  a  massive 
Galactic  halo. 

In  general,  all  the  simulations  for  the  Magellanic  Clouds  (Lin 
&Lynden-Bell  1982;  Gardiner  etal.  1994;  Moore  &  Davis  1994; 
Heller  &  Rohlfs  1994;  Lin  et  al.  1995;  Gardiner  &  Noguchi 
1996;  Yoshizawa  &  Noguchi  2003;  Mastropietro  et  al.  2005; 
Bekki  &  Chiba  2005;  Connors  et  al.  2006;  Besla  et  al.  2007) 
are  based  on  the  backwards  integration  of  the  equations  of 
motion,  first  applied  by  Murai  &  Fujimoto  (1980).  Most  papers 
consider  an  isothermal  sphere  with  a  given  constant  rotational 
velocity  at  a  large  galactocentric  distance,  while  three  of  the 
most  recent  calculations  use  NFW  dark  matter  halo  profiles. 
Not  surprisingly,  the  mass  (profile  and  amount)  of  the  Galaxy 
is  an  important  source  of  uncertainty  in  the  orbital  models.  The 
masses  of  the  Clouds  are  as  well  another  source  of  error,  their 
current  distorted  state  makes  any  dynamical  estimate  of  the  mass 
a  difficult  task. 

Given  the  large  mass  of  the  Clouds,  dynamical  friction15 
can  significantly  reduce  their  perigalacticon  distance  as  they 
move  through  the  Galactic  halo.  A-body  simulations  of  self- 
gravitating  particles  naturally  consider  it  by  default,  while  in 
other  cases,  it  is  accounted  for  by  using  an  analytical  expression 
given  by  Binney  &  Tremaine  (1987).  Recent  studies  (Hashimoto 
et  al.  2003;  Just  &  Penarrubia  2005)  have  found  though,  that 
the  latter  tends  to  circularize  orbits  to  excess  when  compared 
to  equivalent  A-body  simulations,  thus  some  simulations  have 
scaled  down  its  effects. 

The  important  point  to  consider  here  is  that  a  high  orbital 
eccentricity,  or  equivalently  a  high  transverse  motion  with 
respect  to  the  Galaxy,  is  needed  to  lead  to  the  formation 
of  the  high-velocity  Magellanic  Stream.  Therefore,  the  initial 
conditions  of  the  Clouds’  orbits  must  be  such  that  even  after 
the  effects  of  dynamical  friction,  orbital  eccentricity  is  high 
enough  for  the  Stream  to  be  formed.  Interestingly  enough,  Lin 
et  al.  (1995)  had  pointed  out  that  a  hyperbolic  encounter  of  the 
Clouds  with  the  Galaxy,  in  which  they  are  passing  by  for  the 
first  (and  only)  time,  could  lead  to  the  tidal  stripping  of  gas 
segments  that  would  later  infall  and  trail  rapidly  behind.  Such 
a  model  was  discarded  then,  on  the  basis  that  the  Clouds  were 
presumed  bound  to  the  Galaxy. 

Another  source  of  dynamical  friction,  the  LMC  dark  matter 
halo,  is  considered  as  well  in  the  calculations  of  the  Magellanic 
orbits  (Bekki  &  Chiba  2005).  Its  effects  are  quite  important  on 
the  binarity  of  the  galaxy  pair,  exerting  a  significant  frictional 
drag  on  the  SMC  when  it  penetrates  the  LMC  halo.  If  they 

14  By  specifying  the  orbital  inclination  of  the  LMC  with  respect  to  the 
Galactic  plane  and  its  perigalactic  distance,  the  equations  of  energy  and 
angular  momentum  conservation  may  be  solved  to  yield  the  three  components 
of  the  LMC  space  velocity. 

15  The  “retardation”  of  a  moving  object  when  it  passes  through  a  region  with 
non-vanishing  mass  density,  caused  by  its  gravitational  interaction  with  the 
particles  of  that  continuous  mass  (Chandrasekhar  1943). 
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get  close  enough,  they  would  merge  quickly,  and  since  this  has 
not  occurred  yet,  it  implies  either  that  this  force  is  negligible 
or  that  the  Clouds  became  bound  to  each  other  relatively 
recently  (Gardiner  et  al.  1994).  Bekki  &  Chiba  (2005)  also 
found  that  under  these  conditions,  the  Clouds  cannot  keep  their 
binary  status  for  more  than  Gyr  in  the  past.  This  opens  the 
possibility  of  the  Clouds  being  coupled  only  recently. 

Curiously,  the  very  first  dynamical  model  of  the  Magellanic 
System  (Murai  &  Fujimoto  1980)  had  to  make  an  extensive 
search  for  the  binary  state  that  could  produce  the  Magellanic 
Stream.  The  choice  of  tangential  velocities  for  SMC  had  to 
be  so  specific,  that  they  could — theoretically — predict  it  within 
±5  km  s'1.  This  prompted  them  and  others  to  make  the  first 
inferred  estimates  of  the  proper  motion  of  the  Clouds.  This  early 
model  though,  did  not  include  a  dark  matter  halo  for  the  LMC. 

A  similar  situation  was  faced  later  on  by  Heller  &  Rohlfs 
(1994),  who  explain  that  if  the  perturbing  forces  acting  on  the 
LMC  and/or  SMC  have  to  be  small  enough  to  leave  the  binary 
system  intact  while  simultaneously  producing  long  streams 
similar  in  shape  to  the  Stream,  then  the  evolution  time  of  these 
streams  has  to  be  rather  long  and  very  special,  and  properly 
chosen  initial  conditions  for  the  test  particles  in  the  simulations 
have  to  be  adopted. 

One  point  in  which  all  orbital  models  agree,  is  that  the  Clouds 
had  a  recent  encounter,  sometime  between  200  to  500  Myr 
ago.  In  fact,  several  models  also  found  that  the  binarity  of  the 
Clouds,  regardless  of  how  long  it  has  been  in  place,  was  most 
probably  broken  at  this  last  collision,  which  happened  very  close 
to  their  perigalacticton  distance.  In  other  words,  the  Milky  Way’s 
powerful  gravitational  tide  has  disrupted  the  pair. 

In  general,  searching  for  the  appropriate  binary-bound  orbits 
for  the  Clouds  that  can  later  be  used  to  reproduce  the  Magellanic 
System,  has  been  a  difficult  fine  tuning  task.  The  challenge  is 
now  even  harder,  since  the  most  recent  HST  measurements  of 
the  LMC  and  SMC  proper  motion  (Kallivayalil  et  al.  2006b, 
2006a;  Piatek  et  al.  2008)  seem  to  indicate  that  the  Clouds  are 
traveling  too  fast  to  have  ever  been  bound  to  the  Milky  Way. 

Numerical  simulations  of  the  LMC  orbit  by  Besla  et  al. 
(2007),  based  on  these  new  numbers,  suggest  that  the  Large 
Cloud  is  “plunging”  in  a  highly  eccentric  parabolic  orbit,  on  its 
first  passage  about  the  Milky  Way.  At  such  speed,  dynamical 
friction  is  negligible,  but  the  choice  of  the  Galactic  potential 
(isothermal  sphere  versus  NFW)  introduces  dramatic  changes 
in  the  orbital  history  of  the  LMC.  In  an  isothermal  sphere, 
the  LMC  has  indeed  a  bound  orbit,  although  with  an  increased 
period  and  apogalacticon  distance  compared  to  previous  models. 
In  an  NFW  profile,  the  “best  case  scenario”  of  a  bound  orbit  has 
a  period  of  about  a  Hubble  time,  and  reaches  an  apogalacticon 
distance  of  ^550  kpc. 

Despite  such  three-dimensional  differences,  their  projected 
orbits  on  the  sky  were  in  good  agreement,  so  their  predicted 
location  compared  to  the  Magellanic  Stream’s  great  circle,  could 
not  be  used  to  distinguish  Milky  Way  mass  profile  models.  More 
importantly  though,  in  both  cases  the  projected  orbit  did  not 
trace  the  Magellanic  Stream,  deviating  from  it  by  about  10°  on 
the  sky.  Adding  the  SMC  into  the  calculations  did  not  reduce 
the  disagreement,  nor  using  the  weighted  average  of  all  pre- 
HST  proper  motions.  Besla  et  al.  (2007)  argue  then  that  the 
usual  criteria  to  validate  a  dynamical  model  of  the  Magellanic 
System,  its  ability  to  match  an  orbit  with  the  Stream,  is  no  longer 
acceptable. 

The  conclusion  of  Besla  et  al.  (2007),  that  the  Clouds  are 
on  their  first  passage  relies  heavily  on  the  large  value  for 


jLa  cos  5  that  was  measured  with  HST  for  the  LMC.  Lin  et  al. 
(1995)  and  others  have  in  fact  argued  that  it  is  hard  to  explain 
how  a  bound  LMC  could  have  maintained  a  high  angular 
momentum  perpendicular  to  the  Galaxy’s  rotation  axis  for  so 
long.  Raychaudhury  &  Lynden-Bell  (1989),  Shuter  (1992)  and 
Byrd  et  al.  (1994)  have  suggested  that  an  early  interaction  with 
M31  could  be  the  source  of  such  a  high  tangential  velocity.  In 
any  case,  such  a  condition  is  easier  to  understand  with  an  LMC 
not  bound  to  the  Milky  Way.  In  addition,  the  Clouds  are  the  only 
gas-rich  dwarf  galaxies  at  small  galactrocentric  distances  (van 
den  Bergh  2006),  different  from  the  rest  in  the  Local  Group.  All 
these  reasons  are  used  by  Besla  et  al.  (2007)  to  support  the  case 
of  an  unbound  LMC. 

Notably,  all  LMC  proper  motion  results  after  Kallivayalil 
et  al.  (2006b),  including  this  work,  have  produced  lower  values 
of  fAa  cos  8.  Therefore,  a  bound  orbit  with  a  long  period  and  a 
large  apogalacticon  distance  is  still  a  scenario  compatible  with 
the  most  recent  results.  As  for  the  binarity  of  the  Clouds,  a  recent 
period  of  joint  orbits  could  be  enough  to  explain  the  common 
features  between  the  Clouds,  that  is  the  H  i  envelope,  the  Bridge 
and  even  the  Stream,  and  can  account  for  the  different  star 
formation  history  and  chemical  evolution  in  each  Cloud,  that 
point  to  a  separate  origin  and  place  of  birth. 

To  conclude,  the  search  for  a  realistic  orbit  of  the  Magellanic 
Clouds  is  far  from  over.  The  space  velocities  obtained  in 
this  investigation  are  supportive  of  a  scenario  in  which  the 
Magellanic  Clouds  are  possibly  currently  unbound  from  each 
other,  with  the  LMC  traveling  at  a  velocity  that  is  high  enough 
to  make  it  nearly  unbound  to  the  Galaxy.  But  having  nowadays 
a  much  more  precise  measurement  of  the  proper  motion  of  the 
Clouds,  has  not  facilitated  our  understanding  of  their  dynamics 
and  has  instead  opened  new  questions  and  placed  all  the 
constructed  models  in  doubt. 

9.  CONCLUSIONS  AND  FUTURE  WORK 

A  catalog  of  absolute  proper  motions  containing  1,448,438 
objects  has  been  obtained  from  SPM  material,  supplemented 
with  UCAC2  data.  The  catalog  covers  an  estimated  area  of 
450  deg2  except  for  the  inner  regions  of  LMC,  SMC,  and  47 
Tucanae,  where  the  high  stellar  density  made  it  impossible  to 
obtain  accurate  cross-identification  of  the  stars. 

Samples  of  3822  LMC  stars  and  964  stars  were  selected  from 
the  catalog  to  measure  the  mean  proper  motion  of  the  Magellanic 
Clouds.  The  results  obtained  are 

(fia  cos  5,  /t.s)lmc  =  (+L89,  +0.39)  ±  (0.27, 0.27)  mas  yr-1 
and 

(tia  cos  <5,  /rj)sMC  —  (+0.98,  —1.01)  ±  (0.30, 0.29)  mas  yr-1 . 

Our  much  more  precise  relative  proper  motions  with  respect  to 
the  photometrically  selected  Galactic  Disk  dwarf  stars,  enabled 
us  to  obtain  the  proper  motion  of  the  SMC  with  respect  to  the 
LMC,  with  significantly  smaller  uncertainties: 

(Uccms,  m)smc— LMC  =  (-0.91,  -1.49)  db  (0.16, 0.15)  mas  yr-1. 

This  was  used  to  obtain  new  independent  and  more  precise 
proper  motions  for  the  SMC,  based  on  the  more  accurate  LMC 
proper  motions  of  other  authors.  It  was  also  used  to  confirm  if 
their  separate  measurements  of  the  SMC  proper  motion  were 
consistent  or  not  with  our  results. 

After  a  comparison  in  an  absolute  and  relative  sense  with 
previous  proper  motion  results,  followed  by  a  discussion  of  the 


No.  6, 2010 


PROPER-MOTION  STUDY  OF  THE  MAGELLANIC  CLOUDS  USING  SPM  MATERIAL 


1949 


orbital  models  of  the  Magellanic  Clouds  based  on  those  results, 
we  conclude  that  our  proper  motions  are  compatible  with  the 
LMC  and  SMC  being  bom  and  formed  as  separate  entities, 
which  later  joined  in  a  temporary  binary  state  for  the  past 
few  Gigayears,  being  recently  disrupted  by  the  Milky  Way  in 
their  most  recent  perigalaticon  passage  about  200  Myr  ago.  The 
Clouds  orbits  are  marginally  bound  to  the  Milky  Way,  possibly 
following  a  very  elongated  but  still  periodic  orbit  around  the 
Galaxy. 

The  search  for  a  realistic  orbit  of  the  Magellanic  Clouds 
is  far  from  over.  Having  (formally)  very  accurate  and  precise 
space-based  proper  motions  for  the  Clouds  has  not  facilitated 
our  understanding  of  their  dynamics  but  has,  instead,  opened 
new  questions  and  has  placed  all  dynamic  scenarios  of  the 
Magellanic  System  in  doubt.  As  of  today,  it  is  still  unclear  if 
the  Magellanic  Stream  and  the  Leading  Arm  are  caused  mostly 
by  a  tidal  interaction  or  are  the  result  of  the  ram  pressure  of  the 
Galactic  halo  on  the  gas  of  the  Clouds. 

Given  the  inherent  difficulties  in  measuring  an  accurate  proper 
motion  for  the  Magellanic  Clouds,  the  obvious  dangers  that 
systematic  errors  pose  in  those  measurements  and  the  fact  that 
the  dynamical  models  of  the  Magellanic  System  are  extremely 
sensitive  to  small  variations  in  the  proper  motions  of  the  Clouds, 
we  believe  that  we  are  not  yet  in  the  position  of  considering  them 
known  parameters  in  the  orbital  calculation.  But  we  are  getting 
closer. 

9.1.  Future  Work 

From  the  very  beginning,  this  investigation  of  the  proper 
motion  of  the  Magellanic  Clouds  was  known  to  be  restricted 
by  the  conditions  and  characteristics  of  the  SPM  material.  We 
have  proven  that  even  under  those  constraints,  these  data  are 
able  to  produce  independent  significant  results  on  the  proper 
motion  of  nearby  dwarf  galaxies.  Thus,  there  is  certainly  room 
for  improvement. 

Second  epoch  CCD  data  over  the  whole  field  are  needed, 
especially  in  the  inner  areas  of  the  Clouds.  Since  the  definition 
of  an  adequate  uncontaminated  relative  reference  frame  in  those 
areas  is  difficult,  it  is  also  necessary  to  devise  an  adequate  proper 
motion  reduction  method,  to  precisely  link  these  fields  into  the 
general  global  relative  reference  frame.  A  scaled-down  version 
of  the  block  adjustment  method  could  be  put  in  place  here,  in 
which  plate-size  fields  of  view  are  assembled  first,  and  then 
another  solution  is  run  to  paste  those  into  the  larger  global 
system. 

Since  the  dominant  limiting  factor  in  terms  of  precision 
is  the  plate  measurement  errors,  additional  improvement  can 
be  achieved  by  re-scanning  the  first-epoch  plates  used  in  this 
work  with  the  Yale  PDS,  which  yields  two  times  more  precise 
positions  than  the  currently  measured  ones.  As  PDS  scanning 
of  a  full  plate  is  very  time  consuming,  a  subset  of  properly 
pre-selected  objects  should  be  measured. 

Although  additional  HST  follow  up  observations  are  already 
planned  to  improve  the  space-based  proper  motions,  and  other 
research  teams  are  still  working  on  additional  CCD  ground- 
based  measurements,  these  studies  are  still  limited  since  they 
must  correct  their  observed  proper  motions  into  a  center-of- 
mass  value.  Therefore,  it  is  still  worth  trying  to  improve  our 
results,  since  they  offer  a  wide  field  extended  coverage  of  the 
Magellanic  Clouds. 

Another  future  work  being  considered,  is  to  search  in  the 
intercloud  region  for  coherent  structures  in  the  relative  proper 
motion  space,  to  identify  stars  whose  motion  is  directed  towards 


the  LMC.  This  will  require  us  to  find  what  possibly  is  a 
small  number  of  stars  spread  over  an  extended  area  that  share 
systematic  proper  motions. 
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